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The deflection of mesons in a magnetized ferromagnetic medium was investigated. A beam 
of mesons was made to pass through 9 cm of iron, and the resulting distribution of the beam 
was observed. Two arrangements were employed. In the first arrangement, the deflection due 
to the field caused a fraction of the mesons to hit a counter placed out of line with the others. 
An increase of sixty percent in the number of coincidences was recorded when the iron was 
magnetized. In the second arrangement, all the counters were arranged in line, and the de- 
flection due to the field caused an eight percent decrease in the number of coincidences. These 
results are compared with theoretical predictions deduced from the known momentum spectrum 
of the mesons and from the geometry of the arrangement. The observed effects agree as well 
as can be expected with those calculated under the assumptions that the effective vector inside 
the ferromagnetic medium is the induction B, and that the number of low energy mesons is 
correctly given by the range-momentum relation. 





INTRODUCTION 


HE idea of observing the deflection of a 

cosmic-ray particle passing through a mag- 
netized ferromagnetic medium was independently 
suggested and carried out experimentally by 
Rossi! and Mott-Smith.? Both these experi- 
menters obtained negative results, which were 
interpreted at the time as meaning that the 
energy of the cosmic-ray particles was higher 
than expected. 

These experiments provoked a discussion 
among theoretical physicists about the effective 
vector acting on a charged particle moving 
through a ferromagnetic medium. The classical 
theory gives a result that depends on the par- 
ticular model chosen for the elementary magnets.* 


(1930) Rossi, Nature 128, 300 (1931); Lincei Rend. 2, 478 
on LM. Mott-Smith, Phys. Rev. 37, 1001 (1931) ; 39, 403 
3 W. F. G. Swann, Phys. Rev. 49, 574 (1936). 


Moreover, one can hardly rely on the results of 
the classical theory since ferromagnetism is due 
to the electron spin and hence is a typically non- 
classical effect. Fortunately, a calculation based 
on Dirac’s theory of the electron has been per- 
formed by Weizsaecker,‘ the unambiguous result 
being that the effective vector acting on the 
particle is the magnetic induction B. 

Further experiments were performed by AIl- 
varez,® who investigated the deflections of 8-par- 
ticles passing through a thin magnetized iron 
sheet. The observed effect was much smaller 
than expected, and he concluded that the effec- 
tive vector is at most of the order of one-tenth 
of B. ' 

Danforth and Swann® performed another ex- 
periment on cosmic rays. Their counters were 
*C. F. v. Weizsaecker, Ann. d. Physik 17, 869 (1933). 

5L. Alvarez, Phys. Rev. 45, 225 (1934). 


®*W. E. Danforth and W. F. G. Swann, Phys. Rev. 49, 
582 (41936). 
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Fic. 1. Counter arrangement. Counter D is represented 
in the position outside of the beam defined by counters B 
and C. Cross-marked counters are used as anticoincidence 
counters. 


arranged in a straight line, an iron block being 
placed between the two last counters. A decrease 
of the order of ten percent in the number of 
coincidences was observed when the iron was 
magnetized. This result was interpreted, on the 
basis of the then current ideas about the nature 
and energy distribution of the cosmic-ray par- 
ticles, as meaning that the effective field was 
of the order of B. 

The preceding summary of the results of earlier 
experimenters shows that it was worth while to 
undertake a new investigation of the problem 
in view of the improved technique now avail- 
able. Moreover, the present knowledge of the 
nature of cosmic-ray particles at sea level affords 
an unambiguous interpretation of the results. 


EXPERIMENTAL ARRANGEMENT AND RESULTS 


In order to obtain results which admit of a 
clear interpretation, only one kind of particle 
must be selected. Moreover, radiation processes 
and the consequent cascade multiplication make 
the soft component wholly unfit for experiments 
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of this sort. Therefore, one must employ a beam 
of mesons selected by preliminary absorption of 
the soft component. 

The choice of the optimum thickness of mag- 
netized iron to be employed requires some con- 
sideration. For a given thickness /, there is a 
minimum value of the momentum Pmin of the 
mesons that are able to traverse the iron. For a 
given momentum ?, and hence for a given radius 
of curvature R of the path, the deflection wil] 
be approximately proportional to /, and in- 
versely proportional to p. Were the momentum 
proportional to the range, the maximum angular 
deflection obtainable would be independent of 
the thickness. However, in the low energy region 
the range increases rapidly with momentum, 
and it is only in the region of highly relativistic 
energies that the range-momentum relation be- 
comes linear. Consequently, the thickness of 
iron must correspond to particles whose kinetic 
energy is at least as large as the self-energy. On 
the other hand, the use of an excessive thickness 
of iron would excessively reduce the intensity of 
the beam, owing to geometrical factors. When 
all these factors are considered, it is seen that 
the most favorable thickness of iron is of the 
order of ten cm. In the actual experiments, a 
thickness of 9 cm was employed, corresponding’ 
to a minimum momentum of 2.1X10* ev/c 
(throughout this paper, the unit of momentum 
will be 108 ev/c). 

Due allowance must be made for the scattering 
of mesons by the iron. Since the average angle 
for multiple Coulomb scattering and the mag- 
netic deflection are, for relativistic energies, both 
inversely proportional to the energy, scattering 
will equally affect the results at all energies, 
provided these are sufficiently high. Hence the 
effect of scattering, although important, does 
not appreciably affect the choice of the optimum 
thickness. 

The arrangement of the counters is shown in 
Fig. 1. Counters A, B, and C, with the interposed 
ten cm of lead, select a beam of mesons. A fourth 
counter D can be placed either in line with the 
other three, or wholly outside of the beam defined 
by the coincidences (BC). The iron block can 
be placed between the counters C and D. Two 


7B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
(1941). 
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anticoincidence counters E are placed one on 
each side of the counter D, the purpose being to 
discriminate against showers that might give 
fourfold coincidences (ABCD) and against par- 
ticles that, traversing the iron, produce second- 
aries. By displacing the counter D, the definition 
of the meson beam was tested. It was found 
that when the counter D was just outside of 
the beam, the number of fourfold coincidences 
dropped to about one percent of the number 
recorded with the four counters in line. This 
very low “background” is mainly responsible for 
the success of that part of the experiment in 
which the counters were placed out of line and 
was made possible by the use of anticoincidence 
counters. 

With counter D out of line, the number of 
fourfold coincidences was almost doubled when 
the iron block was interposed, owing to scattering 
by the iron. This effect, however, did not prevent 
the measurement of the deflection due to the 
magnetic field, but made it necessary to extend 
the measurements over a longer period. The 
readings were alternately taken with and without 
field, and the difference was assumed to be due to 
the effect of the field. 

The iron block, 94.5 cm? in cross section and 
25 cm in length, constituted part of a closed 
magnetic circuit and was kept magnetized at 
15,000 gauss by an electromagnet dissipating 
40 watts. 

Counters A, D, and E had an effective diam- 
eter of 2.9 cm and an effective length of 37 cm. 
Counters B and C had an effective diameter of 
2.3 cm and an effective length of 22 cm. The 


TasBLeE I. Fourfold coincidences with magnetized and 
unmagnetized iron. Counters out of line. 








Coincidences 
per hour 


Number of 
coincidences 


Number of 


Magnetic 
hours 


induction B 





0.99+0.03 
1.65+0.05 
1.58+0.05 


0 1200 
+ 15,000 650 
— 15,000 650 


1186 
1073 
1028 








distances between the counters are shown to 
scale in Fig. 1. 

The results obtained with the counters out of 
line are summarized in Table I. All the counts 
were taken with counter D placed 0.4 cm out- 
side of the beam defined by the threefold coin- 


cidences (ABC). It is seen that the magnetic 
field increases the number of counts by about 
sixty percent, a far larger effect than any pre- 
vious experimenter has been able to obtain. The 
sign of the field is indicated as positive when the 


TaBLE II. Fourfold coincidences with magnetized and 
unmagnetized iron. Counters in line. 








Magnetic Number 
induction Number of coin- 
of hours cidences 


0 450 13321 
15,000 450 12368 


Coincidences Difference Difference 
per hour - per hour percent 


29.6+0.26 
27.540.25 2.140.36 7.141.2 











positively charged particles are deflected into 
counter D. The difference between the numbers 
of positive and negative particles is within the 
statistical error. There is, however, a slight 
indication of an excess of positive mesons, as 
reported in the literature.*® 

The experiment with the counters in line 
simply consisted in counting the fourfold coin- 
cidences with and without magnetic field. Ob- 
viously, in this case there is no point in reversing 
the field, as the arrangement is symmetrical with 
respect to positive and negative mesons. The 
decrease in the number of fourfold coincidences 
due to the field was of the same order of magni- 
tude as that observed by Danforth and Swann. 
The results are summarized in Table II. 


INTERPRETATION OF RESULTS 


We shall first discuss the results obtained with 
the counters out of line. 

In order to compare the experimental results 
with the theory, one must know the differential 
momentum spectrum of the mesons f(p) and the 
probability w(p) that a meson recorded by 
counters B and C will be deflected into counter D 
by the action of the field, assuming that the 
induction B is effective. Let the function f(p) 
be normalized so that 


J “fedp=1, 


then the ratio of the number of fourfold coin- 
cidences due to the field to the number of incident 
mesons [coincidences (BC)] will obviously be 


® H. Jones, Rev. Mod. Phys. 11, 235 (1939). 














given by: 
n= {  w(p)fiprap. 


min 


The spectral distribution f(p) can be taken 
from data available in the literature, and the 
probability w(p) can be calculated from the 
geometry of the experiment. To perform this 
calculation, we shall make the following simplifi- 
cations: (a) We shall assume that the meson 
spectrum is not appreciably affected by filtering 
through the roof and two floors of the building, 
as well as through the lead block. The distortion 
in the energy distribution thus introduced is 
certainly small, and probably within the limits 
of accuracy of the present knowledge of the 
spectrum. (b) The deflections due to scattering 
and to the magnetic field in the iron block are 
assumed to be independent of each other; that 
is, the effect of the field is obtained by sub- 
tracting the background observed with the un- 
magnetized iron. Actually, scattering will some- 
what alter the intensity distribution in the 
undeflected beam, especially for low momentum 
mesons that are mainly responsible for the 
magnetic effect. This is believed to represent the 
main source of error in comparing the experi- 
mental results with the theory. (c) The radius of 
curvature of the meson path varies as the meson 
loses energy in the iron. To simplify the calcula- 
tions, an average radius will be substituted. 

The momentum spectrum of the mesons is 
wholly determined, for low energies, by the range- 
momentum relation. This relation, if we assume 
that the whole energy loss is due to collisions, 
is well known theoretically. The resulting mo- 
mentum spectrum has been calculated and repre- 
sented graphically, with an approximation suffi- 
cient for our purposes, by Euler and Heisenberg.® 
The spectrum at high energies has been deduced 
by several experimenters from the measurement 
of the curvature of tracks in the cloud chamber. 
A composite spectrum was obtained by joining, 
at p=14, the theoretical curve of Euler and 
Heisenberg with the experimental data of Jones.*® 
This procedure was necessary in order to normal- 
ize the function f(p) although mesons with 
momentum larger than 14 are ineffective in the 


* H. Euler and W. Heisenberg, Ergeb. d. exakt. Naturwiss. 
17, 1 (1938). 
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present experiment. The differential spectrum 
used is shown graphically in Fig. 2. 

The probability w(p) for a meson recorded by 
counters B and C to be deflected into counter D 
was calculated as follows. The intensity distri. 
bution of the undeflected beam in the plane of 
the counter D is readily evaluated geometrically, 
To a first approximation, the effect of the field 
on the particles of one sign can be taken as a 
shift of the whole beam toward one side, by an 
amount inversely proportional to the radius of 
curvature and hence to the momentum.. The 
fraction of the shifted beam that passes through 
counter D is calculated. 

The function w(p) is evaluated for p between 
0 and 14 (for higher values of p the probability 
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Fic. 2. f(p), differential momentum spectrum; multiply 
ordinates by 0.1. w(p), probability for a meson to be 
deflected into counter D. f(p)w(p), product of the two 
functions; multiply ordinates by 0.02. 







is practically zero) and is represented in Fig. 2, 
The figure also represents the product f(p)w(p); 
from this latter curve it is seen that most of the 
contribution to the magnetic effect comes from 
mesons that have little more than the energy 
necessary to traverse the iron block. 

A numerical integration gives: 











14 


_ f w(p)f(p)dp = 0.028. 






In order to compare this figure with the experi- 
mental data, it is also necessary to know what 
fraction N of the beam defined by the coin- 
cidences (BC) passes through counter D when 
the latter is in line with the other two. This 
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fraction can immediately be evaluated from the 
geometry and is 0.49. The theoretical value of 
the ratio 2/N is therefore: 


(n/N) theor = 0.057. 


Since only mesons of oné sign can be deflected 
into counter D, we shall take the sum of the 
effects due to the positive and negative fields 
after subtracting the effect measured without 
field. This sum is 1.25+0.08 coincidences per 
hour. The number of fourfold coincidences with 
the counters in line is 29.6 per hour. Hence 


(n/N exper = 1.25/29.6 = 0.042+-0.003. 


Considering the uncertain factors involved in 
the calculations, we may state that there is no 
significant discrepancy between the observed 
and the calculated effect. 

We shall now consider the experiment with the 
counters in line. The calculations are similar to 
those discussed in the preceding case; hence the 
procedure followed will be briefly outlined. 

In this case, the displacement of the beam due 
to the field gives rise to a decrease in the number 
of particles recorded by counter D. This decrease 
was calculated, from the geometry of the beam, 
as a function of momentum, and a numerical 
integral was again evaluated on the basis of the 
same spectral distribution as before. The theo- 
retical result is that the field should produce a 
decrease of 8 percent in the number of fourfold 
coincidences. From Table II it is seen that the 
observed effect is 7.1+1.2 percent. Hence, the 
experimental and theoretical values agree to 
within the statistical error. 

We now wish to discuss the significance of the 
present experiments. The agreement between the 
calculated and the observed effects was obtained 
under two assumptions: (a) the effective vector 
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acting on the particles inside the magnetized 
medium is the induction B; (b) the number of 
low energy mesons is correctly given by the 
range-momentum relation. Were the effective 
field considerably smaller than B, or were the 
number of low energy mesons smaller than 
assumed, then the observed effect would not 
have agreed with the theoretical predictions. 

It would at first appear that the present experi- 
ment could not be considered as a test of both of 
the above-mentioned assumptions. We must con- 
sider, however, that if either of the two deter- 
mining factors disagreed with the theoretical 
predictions, this disagreement would be only in 
one direction. The effective vector acting on a 
particle could conceivably be smaller than B 
(under certain classical assumptions, for ex- 
ample, one obtains values that are intermediate 
between H and B); but one could hardly imagine 
a reasonable hypothesis under which the effec- 
tive vector could be larger than B. Regarding 
the number of low energy mesons, here again 
one may readily imagine mechanisms!® capable 
of removing mesons from the beam, and hence 
of rendering the number of slow mesons smaller 
than assumed; but it would be difficult to con- 
ceive of an effect that could make the number of 
slow mesons higher than calculated from the 
range-momentum relation. We conclude that the 
agreement of the experimental results with the 
theory is a good proof, at least, of the approxi- 
mate validity of both hypotheses, namely, that 
the effective vector is B and that the number of 
slow mesons is correctly given by the range- 
momentum relation. 

10 An interaction that might produce this effect has been 
discussed by J. Hamilton, W. Heitler, and H. W. Peng, 
Phys. Rev. 64, 78 (1943). Their theoretical conclusions do 


not seem to be supported by the results of the present 
experiments. 
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PHYSICAL REVIEW 


NUMBERS 1 AND 2 


The Effects of Various Ions on the Infra-Red Absorption of Water 


DubLey WILLIAMS, University of Oklahoma, Norman, Oklahoma 
AND 
WALTER MILLETT, University of Florida, Gainesville, Florida 
(Received May 1, 1944) 


The effects of various non-hydrolyzing salts on the infra-red spectrum of water have been 
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studied in the region between 1 and 7,4. Variations in the positions of the absorption maxima at 
1.46y, 1.79, 1.96u, 3.04, 4.74, and 6.14 have been observed. The changes produced by the 
solutes have been compared with existing data on the effects of temperature variation on 
the absorption of water. The results indicate that, in general, large ions produce changes 
similar to the changes produced by increasing the temperature of the pure solvent, whereas 
small ions produce the reverse effect. These results are in agreement with data on the other 


physical properties of water and aqueous solutions. 








N the basis of x-ray diffraction data, Bernal 
and Fowler! have developed a theory of the 
intermolecular structure of water and aqueous 
solutions. According to this theory, the molecules 
of water in the liquid state are arranged in local 
crystalline lattices. The extent of these lattices 
and the regularity of the arrangement of the 
molecules within the lattices are greatest at low 
temperatures. The effects of the presence of ions 
in the water depend upon the polarizing powers 
of the ions, polarizing power being defined as 
the ratio of ionic charge to ionic radius. In a 
sequence of ions of equal charge, the larger ions 
decrease the regularity of molecular arrange- 
ment and produce changes similar to those pro- 
duced by elevation in the temperature of pure 
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Fic. 1. The transmission spectrum of water in the vicinity 
of the 1.79-» band. 


1J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 
(1933). 


water ; smaller ions tend to enter into the crystal- 
line lattices and preserve regularity. 

As the infra-red absorption bands of water 
have been found to change in position as the 
temperature is varied,** it has been of interest 
to study the influence of various ions on these 
bands in order to compare the observed effects 
with the effects predicted on the basis of the 
Bernal-Fowler theory. Such studies have already 
been made for the 4.74 band.** The purpose of 
the present investigation has been to extend 
these studies to include all the major infra-red 
bands between iy and 7yz. 

The spectrometer used was an instrument of 
the Wadsworth-Littrow type equipped with a 
large glass prism for use between 1y and 2y and 
a rocksalt prism for use at longer wave-lengths. 
Absorption cells with glass windows were used 
at wave-lengths shorter than 2y and fluorite cells 
were used between 2u and 7y. Cell thicknesses 
were adjusted to give optimum measurable ab- 
sorption for each band. The bands studied were 
located at 1.46u, 1.79u, 1.96u, 3.0u, 4.7, and 
6.1u. Assignment of these bands to various types 
of molecular vibration have been made by earlier 
investigators. ® 

The compounds chosen for study were the 
alkali and alkaline-earth halides.’ The ions pro- 

2]. R. Collins, Phys. Rev. 26, 771 (1925). 

*E. Ganz, Ann. d. Physik 28, 445 (1937). 


‘E. K. Plyler and E. S. Barr, J. Chem. Phys. 6, 316 


(1938). 

(1982). Williams and T. N. Gautier, Phys. Rev. 56, 616 
193 

* J. W. Ellis, Phys. Rev. 38, 693 (1931). 

7 All compounds were obtained from Eimer and Amend; 
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SHIFT OF ABSORPTION MAXIMA IN CM” 


Fic. 2. Changes in the characteristic vibrational frequencies of water: The plots show frequency changes produced by 
various positive ions in the presence of common negative ions. 


duced are monatomic and hence have no absorp- 
tion in the region being investigated. In order to 
determine the effects of various cations on a 
given band, a sequence of solutions of compounds 
containing these cations in combination with a 
common anion was studied. The reverse pro- 
cedure was used to ascertain the effects of various 
anions. The same molar concentration was used 
for each solution in a given sequence. In general, 
the concentration was made as large as possible, 
but it was found desirable to use less concentrated 
solutions in the thin cells employed for the 3.0-u 
and 6.1-u bands. 

Definitely measurable shifts produced by 
solutes were observed for all bands except the 
weak band at 1.79». The absorption of water 
in this region is shown in Fig. 1. It will be 
observed from this figure that the 1.79-u band 
appears as a small depression, on the side of 
the intense band with maximum at 1.96u. As 
indicated by Collins,® the location of the 1.79-z 


only small quantities of the cesium and rubidium com- 
pounds were available and these were used for work in the 
short wave-length region. 


8J.R Collins, Fa. Rev. 52, 88 (1937). 


maximum can be determined satisfactorily by 
taking the ratio of the observed transmission to 
an assumed background such as that given by 
the dotted curve in Fig. 1. However, the back- 
ground absorption changes when solutes are 
introduced, since the solutes produce changes in 
the intense 1.96-u band. Hence, it is necessary to 
assume a different background for every solution 
studied and it is felt that the shifts detected by 
such methods are not sufficiently quantitative to 
warrant publication. However, it can be asserted 
that the changes produced by the fluoride and 
iodide ions are of the type to be expected from 
the Bernal-Fowler theory. No definite statements 
can be made about the effects of ions of inter- 
mediate size. 

The results obtained for the other bands are 


. summarized in Figs. 2 and 3. In these figures the 


shifts of absorption maxima are plotted. The 
observed shifts are indicated by the vertices of 
the shaded triangles and the uncertainties in the 
location of the maxima are indicated by the 
length of the bases of the triangles. Concentra- 
tions are given for each sequence of solutions. 
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Fic. 3. Changes in the characteristic vibration frequencies of water: The plots show frequency changes produced by 
various negative ions in the presence of common positive ions. 





The data for water at various temperatures are 
taken from Ganz’ work on the 6.1-y, 4.7-u, and 
3.0-4 bands and Collins’ work on the 1.96-u and 
1.46-u bands. 

Figure 2 gives the shifts produced by various 
positive ions in the presence of common negative 
ions. According to Bernal and Fowler, increasing 
ion size should produce effects similar to those 
produced by increasing the temperature. Of the 
20 sequences summarized in Fig. 2, 14 sequences 
exhibit shifts of the type predicted from the 
theory, 3 sequences show shifts in the opposite 
direction, and 3 sequences show no definite trend. 

In Fig. 3 are shown the shifts produced by 
various negative ions in the presence of common 
positive ions. Of the sequences shown, 16 show 
shifts of the predicted type, one has a shift to 
the opposite direction, and three show no definite 
shift. It should be noted that MgCl. and CaCl, 
produce a definite splitting of the 4.7-u band. 
Six of the sequences for which the shifts were 








not in agreement with those predicted from the 
theory were for data obtained in the 3-y region. 
The broad band in this region is produced by 
three overlapping bands. Hence, it is probably 
desirable to repeat the measurements in this 
region with an instrument of higher resolving 
power and with thinner absorbing layers than 
could be obtained with the type of cells used in 
this study. It would also be desirable to make a 
thorough study of the effects of concentration on 
the changes produced by the solutes. 

The results obtained are in general agreement 
with the predictions made on the basis of the 
Bernal-Fowler theory. However, the fact that 
there are some exceptions probably indicates 
that there are factors other than the polarizing 
power of the ions which should be considered in 
developing a more complete theory. 

The writers wish to express their appreciation 
to Miss Annette Herald for preparation of the 
figures. 
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The Effect of Adsorbed Gases on Contact Electrification* 






Davip E. DEBEAU 
Department of Physics, University of California, Berkeley, California 
(Received May 25, 1944) 









The magnitude of the electrification produced on contact of quartz and of sodium chloride 
on nickel inatmospheres of air, oxygen, and nitrogen has been measured. Reproducible pressure- 
dependent results were obtained that indicate the contact electrification is a surface phe- 
nomenon and depends greatly on the nature of the surfaces involved. At least two layers of 
adsorbed gas play an important part in the phenomenon of contact electrification. In the 
pressure range from 760 to 30 mm at least one layer of adsorbed gas is removed. In this pressure 
range the surface covered with adsorbed gas makes the major contribution to the charge 
separation. At pressures below 0.1 mm at least a second layer of adsorbed gas is removed. 
In this pressure range the major contribution to the charge separation is made by the surface 
from which the second layer of gas has been removed. A fit of the data to Langmuir adsorption 
isotherms was obtained. Water vapor and perhaps hydrogen act as surface poisons in the 
phenomenon of contact electrification. 























ONTACT or frictional electrification and _ in industry would be of greater importance if the 
contact potential difference are two “‘solid phenomenon were better understood. While a 
to solid’’ contact phenomena that lead to the great deal of work has been done, very little 
electrical charging of the solids involved. Contact satisfactory information has been obtained. 
potentials are observed only when both of the Nearly all the published data are qualitative 
solids are metals. This effect may be readily ex- and often are contradictory and not repro- 
plained in terms of the electron theory of metals. ducible.'~' In fact the very nature of the phe- 
The magnitude of the contact potential may nomenon is obscure. It was suggested by Pro- 
either be measured directly or computed from the fessor Loeb that this phenomenon be further 
work functions of the metals. The results thus investigated and that great care be taken to work 
obtained are consistent and reproducible. with reproducible surfaces in a reproducible 
When one or both of the solids is a non-con- environment with particular reference to the 
ductor, the measurement of contact potential is nature and pressure of the gases in contact with 
not possible or rather meaningless. Non-con- the surfaces. 
ductors do not accept or give up electrons readily. Nearly every man who has worked on the con- 
Such bodies do not have the free electrons which tact potentials of metals has noticed and re- 
are necessary for the phenomenon of contact po- marked on the effect of adsorbed gases on the 
tential difference. Therefore the theory which potentials measured. It is now customary to take 
explains contact potentials cannot explain con- great pains to remove all traces of gas, to prepare 
tact electrification. Nor has there yet been pub- clean surfaces, before making a measurement of 
lished a theory which satisfactorily explains this contact potential. 
phenomenon. Recently the importance of adsorbed gas in 
The charging of two non-conductors by friction @nother contact phenomena, friction, has been 


was the first form of electrical charging discovered. Shown.*’ It therefore seems strange that the 
The phenomenon of charging by impact or con- few men who have studied contact electrification 



















tact is a very widespread one encountered in dust ; 
: ° . 1P, E. Shew and C. S. Jex, Proc. Roy. Soc. London 
storms and in many industrial processes where A118, 97-113 (1928). 


there is friction between non-conductors. The 2A. Fleming, Proc. Roy. Soc. London 53, 51-3 (1941). 
control and applications of contact electrification — Proc. Phys. Soc. London 53, 547-53 
Saree ae ‘P. Boning, Zeits. f. tech. Physik 8, 385-98 (1927). 

_ *Submitted to the faculty of the University of California 5H. Israel, Zeits. f. tech. Physik 9, 289-93 (1928). 

in partial fulfillment of the requirements for the degree of ®R. Schnurman, Proc. Phys. Soc. 53, 538-46 (1941). 
Doctor of Philosophy. 7H. Donandt, Chem. Abs. 36, 3466 (1942). 
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recently have paid little or no attention to the 
gases adsorbed by the solids used. 

I have attempted to study quantitatively the 
charge separation produced on the contact of two 
dissimilar solids, one of them a non-metallic, and 
the effect of gas pressure on this charge separation. 


APPARATUS 


The apparatus consists of two nickel flasks 
mounted coaxially and separated by a cylindrical 
insulator made of quartz. One of the flasks is 
filled with small particles of the non-metallic 
solid. The end plates of the flasks are perforated 
so that equilibrium gas pressure may be quickly 
established. The whole is enclosed in a glass 
envelope connected to a vacuum system through 
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Fic. 1. Cross section of the apparatus. 


a ground glass joint. (See Fig. 1 for a cross 
section of the apparatus.) 

When the envelope is inverted, the non-metallic 
particles fall from one flask to the other, carrying 
with them a charge. The potential built up be- 
tween the flasks is measured with a quadrant 
electrometer. From this potential and from the 
known capacity of the apparatus the charge 
separation is computed. 

The total capacity of the system, flasks, 
electrometer, and parallel capacitor is 750yyf. 


MATERIALS 


The quartz particles used were prepared by 
crushing amorphous quartz in an agate mortar. 
The NaCl was prepared from commercial salt. A 
hot saturated solution of commercial salt was 
filtered and then heated until about one-fourth 
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TABLE I. Physical constants of the particles used. 



























Quartz Quartz Quartz NaCl 
Weight 7.02 g 4.74¢ 7.81 ¢g 8.10 g 
Density 1.55 1.55 1.55 2.17 
Total vol. (cc) of non- 
metallic used 4.52 3.05 5.03 3.73 





Av. vol. (cc) of indi- 
vidual particle 
Number of particles 





3.2X107-§ 3.2X10-* 3.2X10-* 3.2x107 
1.41 X105 0.95 X10® 1.57X10®° 1.16 x10 












TABLE II. Quartz on nickel with air in the chamber, 
7.0 grams of quartz. 














Electrometer reading 
















Date Pressure} 1 2 3 4 5 6 Av. Volts 
4/26/43 30 26.7 25.2 25.9 25.8 26.2 26.0 53.3 
4/27/43 31 25.7 25.8 26.6 26.3 25.4 26.1} 26.0 53.3 
7/ 2/43 30 25.3 24.6 24.7 24.0 24.1 24.5 494 
5/ 6/43 7 X10-5| 36.6 35.1 36.5 36.1 35.2 36.2| 36.0 104 
5/ 6/43 264 36.7 35.9 368 36.3 36.1 36.4 107 
4/26/43 42 25.4 24.7 24.9 25.1 24.6 25.0} 24.9 50.5 


























of the water had evaporated. The crystals formed 
were recovered, dried, and crushed in the agate 
mortar. 

Particles of nearly uniform size were segregated 
by differential screening. The average particle 
radius was 0.02 cm. The volume of an average 
particle was computed from the diameter by 
assuming that the particles were spherical. An 
approximate figure for the number of particles 
used was then computed (Table I). 

The gases used were air, oxygen, nitrogen, and 
hydrogen. No attempt was made to purify the 
gases other than to remove water vapor, mercury 
vapor, oil particles, and dust particles. 


EXPERIMENTAL PROCEDURE AND RESULTS 


In order that data taken on different runs be 
comparable, it was necessary that the surfaces 
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Fic. 2. Potential-pressure dependence for the entire 
range studied. 
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EFFECT OF ADSORBED GASES 


used be reproducible. I found that baking the 
apparatus for several days at 200°C in vacuum 
produced nickel, quartz, and NaCl surfaces with 
characteristics that were sufficiently reproducible 
for the needs of this experiment. The apparatus 
was given this treatment before each run. 

For both quartz and NaCl a large number of 
measurements of potential were made at various 
pressures ranging from 7 X 10~* to 760 mm of Hg. 
The rate of increase of pressure between measure- 
ments was slow so that the water vapor and dust 
traps could work efficiently. The rate of decrease 
of pressure between measurements was rapid. 

At any single pressure from five to eight meas- 
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Fic. 3. Potential as a function of pressure—the low 
pressure range. 


urements of the potential were made from a 
similar number of inversions of the apparatus. In 
Table II are listed typical sets of measurements 
taken on different days and at different pressures. 
The difference between the highest and lowest 
reading of a set varies between 2 and 6 percent of 
the average of the set. The differences between 
average values taken for the same pressure on 
different days are of the same order of magnitude. 

The sets of measurements were averaged and 
plotted (Figs. 2 to 6). In all, four sets of runs were 
made. The first run was made with 7.0 g of quartz 
particles in the apparatus. In the second run 8.1 g 
of NaCl were used and in the third run 4.7 g of 
quartz was used. The nickel flasks were then 
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Fic. 4. Quartz: dependence of potential on gas in 
the chamber. 





Nac 


AIR 
02 





| l 
20 40 60 80 
PRESSURE IN MM OF MERCURY 


Fic. 5. NaCl: dependence of potential on gas in 
the chamber. 
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Fic. 6. Potential as a function of pressure—the high 
pressure range. 


cleaned in dilute nitric acid and washed in dis- 
tilled water. The fourth run was made with 7.8 g 
of freshly ground quartz. 

In Fig. 2 potential is plotted against the loga- 
rithm of the pressure in order to give an over-all 
picture of the experimental results. The low 
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Fic. 7. Charge separation in terms of electrons per particle 
as a function of pressure. 


pressure data are plotted in Fig. 3. For conve- 
nience the potential is plotted against McLeod 
gauge readings and thus against a quantity pro- 
portional to the square root of the pressure. In 
Figs. 4-6 the potential is plotted as a function 
of pressure in mm of Hg. 

On the whole the quartz and NaCl behave 
quite similarly in the apparatus. But there was 
one very striking difference. After electrification 
the quartz particles were always negative while 
the NaCl particles were always positive. 

From the capacity of the system and the 
measured potential, I computed the total charge 
separation, and the charge separation in terms of 
electrons per particle (Table III and Fig. 7). 
Although I have given the charge separation in 
terms of the electrons per particle, I do not wish 
to intimate that the mechanism of charging is 
that of a transfer of electrons as such. A possible 
mechanism will be discussed later. 


DISCUSSION OF EXPERIMENTAL RESULTS 
Air, Oxygen, and Nitrogen 


On examining Figs. 4 and 5, we see that for the 
middle pressure range there is no significant 
difference in the behavior of oxygen, nitrogen, 
and their mixture air in this experiment. An ex- 
amination of all of the data available shows that 
there is no real difference in their behavior at any 
pressure that I have studied. 

The high pressure portion of the curve (Fig. 6) 
has a shape similar to that of the Langmuir 
adsorption isotherm.® The analytic expression for 


81, Langmuir, J. Am. Chem. Soc. 38, 2267 (1916); 40, 
1361 (1918). 
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the adsorption isotherm is 
x=aP/(1+aP), (1) 


where x represents the fraction of the surface 
covered by the adsorbed gas or rather the fraction 
of completion of an adsorbed layer. P represents 
the pressure of the gas being adsorbed and a is a 
constant depending on the gas and the surface. If 
the charge separation depends on the presence of 
adsorbed gas, we may thus postulate that the 
potential measured will be proportional to the 
fractional completion of the adsorbed layer. That 
is, 

















V=bx (2) 
or 
V=baP/(1+aP). (3) 
TABLE III. Charge separation. 
Pressure Electrons 
Non- in mm of per 
metal mercury Volts E.s.u. particle 
Quartz 0.97 24.8 55.8 1.21 108 
Quartz 47 67.2 151 3.28 X 10 
Quartz 649 122 273 5.95 X 10 
NaCl 0.81 12.1 27.2 0.72 X10 
NaCl 60 36.5 82.2 2.16 X 10 
NaCl 701 60.5 136 3.59 X 10 
On rearranging, Eq. (3) becomes 
1 1 1 
—=—+-. (4) 
V abP b 


Then 1/V plotted against 1/P should be a 
straight line with slope 1/ab and intercept 1/0, 
The reciprocal values of V and P are plotted in 
Figs. 8-11. The results approximate fairly well to 
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Fic. 8. 7.0 g quartz: plot of the reciprocal of the data for 
determination of adsorption isotherm constants. 










Fic 


ma 
nat 
the 
of 1 
Sin 








(1) 


irface 
.ction 
sents 
aisa 
ce. If 
1Ce of 
t the 
O the 
That 


(2) 


(4) 


be a 

1/b. 
ed in 
ell to 





ta for 








4.75 QUARTZ 














10 Q0! 002 03 


Fic. 9. 4.7 g quartz: plot of the reciprocal of the data for 
determination of adsorption isotherm constants. 


a straight line. The constants obtained by fitting 
straight lines to the points are listed in Table IV. 
These constants may be checked by comparing 
computed and observed values of the potential. 
Such a comparison is made in Tables V and Vi. 
The data used for the comparison were selected 
at random from the available data. The differ- 
ences between the observed and the computed 
values of the potential are of the same order of 
magnitude as the differences between individual 
values of the potential for a given pressure. We 
may safely say that Eq. (3) fits the experimental 
data within the uncertainty of the experiment. 
Let us now inquire into the dependence of the 
constants of Eq. (3) on physical factors. The 
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Fic. 10. 7.8 g quartz: plot of the reciprocal of the data for 
determination of adsorption isotherm constants. 





magnitude of the constant a should depend on the 
nature of the gas adsorbed and on the nature of 
the adsorbing surface, and should be independent 
of the nature of the other solid surface involved. 
Since the magnitude of a did not change when 
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NaCl was substituted for quartz (Table IV), it 
seems likely that the gas involved in the contact 
electrification had been adsorbed by the nickel. 
The other possible but less likely conclusion is 
that the adsorption behavior of quartz and NaCl 
toward the gases used is nearly identical. 

The magnitude of the constant b should depend 
on the nature of both solids involved and on the 
number of particles that can make contact. We 
may assume that the charge separation is pro- 
portional to the area of contact and therefore to 
the number of similar particles that can make 
contact. In the last row of Table IV are listed the 
values for b adjusted to 105 particles. The values 
for the three quartz runs are of the same order of 
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Fic. 11. NaCl: plot of the reciprocal of the data for de- 
termination of adsorption isotherm constants. 


magnitude while the value for NaCl is quite 
different. 

The action of NaCl is different from that of 
quartz not only in that the sign of the charge 
produced is different but also in the magnitude of 
the electrification per particle. 

In arriving at Eq. (3), I assumed that the mag- 
nitude of the charge separation depended only on 
the portion of the surface covered with adsorbed 
gas. If the uncovered portion of the surface also 
makes a contribution to the charge separation, 
Eq. (3) must be changed to-read 


aP dP 
V=b——+e( 1- : (5) 
1+ab 1+dP 


If the uncovered portion of the surface does make 
such a contribution in this pressure range, it was 
so small as to have been obscured by the random 
variations in the results obtained. 
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TABLE IV. Adsorption isotherm constants—the high 
pressure range. 














Quartz 
7.02 ¢ 4.74¢ 7.81 ¢g NaCl 

Slope 0.375 0.64 0.37 0.80 
Intercept 0.0078 0.0133 0.0070 0.0158 
ab 2.67 1.56 2.70 1.25 
b 128.0 76.0 143.0 63.3 
a 0.021 0.021 0.019 0.020 
Number of par- 

ticles 1.41 X108 0.95 <X105 1.57 X10 1.16 X105 
6 adjusted to 10° 

particles 90.9 80.0 91.2 54.4 











At pressures below 30 mm the data no longer 
fit Eq. (3), and at pressures below 0.1 mm the 
magnitude of the charge separation increases 
with decrease in pressure (Figs. 2 and 3). The 
potentials computed from Eq. (3) for pressures 
below 30 mm are lower than the observed po- 
tential. This difference is probably the contri- 
bution of the surface uncovered by the removal 
of the layer of adsorbed gas. The rapid increase in 
potential below 0.1 mm cannot be assigned to 
changes in the first layer of gas, for even at 1.0 
mm but 2 percent of the surface is covered by this 
layer. This increase in potential must be caused 
by further changes in the surfaces making contact. 
Since the change in potential is still pressure de- 
pendent, it is reasonable to assume that the 
change involved is the removal of a second layer 
of adsorbed gas. Since the potential increases 
with decrease in pressure, we may assume that 
the surface uncovered by the removal of this 
second layer is making the major contribution to 
this potential. If this be the case, the dependence 
of potential on pressure should have the form of 


Eq. (6), 
V=c/(1+dP). (6) 


For pressures below 0.1 mm a decent fit of the 
data of Eq. (6) was obtained. The constants for 
the two runs with sufficient low pressure data to 
be so analyzed are listed in Table VII. 

There are then at least two layers of adsorbed 
gas involved in the pressure range studied. At 
atmospheric pressure both layers are essentially 
complete and the charge transfer is large. As the 
pressure is decreased, one of these adsorbed 
layers is removed and the amount of charge 
transferred is reduced. The surface uncovered by 
the removal of the first layer of gas is then of lower 
activity than the same surface covered by the 
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gas. The dependence of potential on pressure in 
this range is given by Eq. (3). At a pressure of 1 
mm the first layer has been nearly removed but 
the second layer has not yet been touched. It is 
in this region that the charge transfer is a mini- 
mum. As the pressure is further reduced, the 
second layer of adsorbed gas is removed uncover- 
ing a surface of higher activity than that of the 
same surface covered by the adsorbed gas. The 
dependence of potential on pressure in this range 
is given by Eq. (6). It is probably a coincidence 
that the charge transfer obtained when both of the 
layers are complete is of the same order of 
magnitude as the charge transfer obtained when 
both layers have been removed. F 

We have not sufficient data to assign the 
second layer to one of the surfaces making con- 
tact. It may even be that what appears to be the 
removal of a single layer of adsorbed gas in each 
of the pressure ranges studied may be the nearly 
coincident removal of a layer of adsorbed gas 
from each of the surfaces. 

Judging from the consistency of the results ob- 
tained with air, oxygen, and nitrogen, I believe 
that I have succeeded in reproducing the solid 
surfaces used by consistent treatment. I cannot, 
however, definitely state the exact nature of the 
surfaces used. It is probable that even at the 


TaBLE V. Comparison of computed and experimental 
values of the potential. 7.0 g quartz. 








Computed Caagyees 





¥ 2.67 P 1+0.021 P Vo-Ve 

4 10.7 1.084 9.9 27.2 +17.3 
11 29.4 1.23 23.9 34.7 +10.8 
25 66.8 1.53 43.6 50.6 + 7.0 
30 80.1 1.63 49.3 53.3 + 4.0 
42 112 1.88 59.4 50.5 — 89 
61 163 2.28 71.5 73.0 + 15 
86 229 2.81 81.0 75.0 — 6.0 
101 269 3.12 86.3 85.0 - 13 
128 342 3.79 90.4 94.0 + 3.6 
147 393 4.09 96.5 96.0 — 0.5 
208 555 5.37 103 103 0.0 
251 670 6.27 106 101 — 5.0 
375 993 8.81 112 118 + 6.0 
463 1235 10.71 115 116 + 1.0 
699 1731 14.61 118 122 + 3.0 








lowest pressures used one or both of the surfaces 
were covered with a layer of adsorbed oxygen or 
nitrogen or both. The work of Roberts® indicates 


* J. K. Roberts, Some Problems in Adsorption (Cambridge 
University Press, 1939). 
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that at pressures below 10-° a metal surface can 
be covered by more than one strongly held layer 
of adsorbed gas. 


Hydrogen and Water Vapor 


In an attempt to use hydrogen in the appa- 
ratus with quartz, I ran into many difficulties. | 
was unable to obtain any reproducible data in the 
hydrogen atmosphere. The charge collected was 
first of the same order of magnitude, then much 
lower than in the case of the other gases. I also 
had considerable difficulty in removing the hydro- 
gen. Baking the apparatus at 200°C for three 
days while it was being evacuated by a mercury 
pump removed only some of the surface hydro- 
gen. A run made in air after this treatment gave 
results that were much too low and which de- 
creased with each inversion of the apparatus. It 
was only after the apparatus had been baked at 
200°C and 20 mm of air for six days that repro- 
ducible results were again obtained. 

The effect of hydrogen on the NaCl-nickel 
system was much the same as its effect on the 
quartz-nickel system. 

The apparatus as it comes from the glass 
blower has been in contact with moist air. In this 
condition extremely low results are obtained for 
the charge separation. The charge separation can 
be returned to the higher reproducible values by 
baking out the apparatus at 200°C and a pressure 
of 7X10-° mm of mercury. 

Apparently water vapor and hydrogen act as 
poisons for the phenomena of contact electrifica- 
tion. The poisoning by hydrogen may have been 
a water vapor effect. It is probable that oxygen 
was never completely removed from the appa- 
ratus before the hydrogen runs. Then the proc- 
esses involved in contact electrification may have 
catalyzed the chemical combination of hydrogen 
and oxygen to form a poisoning water layer on 
one of the surfaces. This mechanism would also 
account for the inconsistent results obtained with 
air when there was some residual hydrogen in the 
apparatus. 


Mechanism of Contact Electrification 


We do not have sufficient data to state defi- 
nitely the mechanism of contact electrification of 
quartz or sodium chloride on nickel. The results 


TaBLeE VI. Comparison of > ety and experimental 
values of the potential. 8.1 g NaCl. 
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TABLE VII. Adsorption isotherm constants—the low 
pressure range. 








7.0 grams quartz 8.1 grams NaCl 





105 70 
17 14 








indicate that the adsorbed gas molecules play an 
important part in the mechanism of charge 
transfer. 

A substance like quartz does not accept or give 
up electrons easily, and therefore a mechanism 
involving the transfer of electrons is unlikely. 
But even amorphous quartz has holes in its short 
range order system from which ions of oxygen are 
missing. Some of these holes are on the surface. 
The presence of these holes should make it easier 
for negative ions to transfer from the nickel to the 
quartz. It is possible that the ions of the adsorbed 
gas molecules, oxygen or nitrogen, are transferred 
from the nickel to the quartz on contact and 
carry the charge. 

In the case of sodium chloride this mechanism 
necessitates either the transfer of positive ions 
from the nickel to the sodium chloride or the 
transfer of negative ions from the sodium chloride 
to the nickel. The transfer of positive oxygen ions 
seems an unlikely process. An alternative process 
is the transfer of negative ions of gas molecules 
adsorbed by the sodium chloride or the transfer 
of chlorine ions from the surface of the sodium 
chloride. In either case the awkward point is the 
constancy of the constant a of Eq. (3) in all the 
experiments. If we accept this latter mechanism 
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and also the constancy of a, we must believe that 
nickel and sodium chloride behave in much the 
same way in adsorbing oxygen and nitrogen. 

Further study of the phenomenon is needed 
before even an hypothesis for the mechanism can 
be made. In spite of the conflicts introduced by 
the ionic mechanism, it still seems more likely 
than an electronic mechanism. 
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An attempt is made to explain the phenomenon of condensation of the gas, on the basis of 
the classical mechanical model of a gas. The condensation of the gas is connected with a 
change in the topological structure of the energy surface in the 6N-dimensional phase space. 
This occurs when the total energy is zero. The agreement with experiment is satisfactory. 


1. INTRODUCTION 


HIS paper is based on the assumption of a 
classical model for a gas: It is composed of 
N point masses attracting and repelling each 
other ; the force between two molecules is derived 
from a potential ar~*—br—, where r is the dis- 
tance between them; a, }, , m are positive con- 
stants (which we will consider as known from 
experiment) with »>m. We neglect internal 
degrees of freedom of the molecules. The differ- 
ential equations governing the motions determine 
a family of trajectories in the 6N-dimensional 
phase space. If the total energy of the system is 
fixed, equal to c, the trajectories are restricted to 
a (6N—1)-dimensional manifold M(c). It is to be 
remembered that, since in a gas there is never 
any information as to the initial conditions of the 
individual molecules, this manifold (and the 
volume, see Section 5) is the only thing that is 
physically given. 

The structure of M(c) is analyzed for varying c 
and it is found that, for c=0, M(c) has one type of 
structure, while for c<0 the structure is quite 
different. Roughly speaking, at c=0, M(c) col- 
lapses from an infinite manifold (c=0) to a finite 
manifold (c<0). It would seem appropriate to 


identify the energy value zero (where the struc. 
ture of the manifold changes) with a state of 
transition of the whole system. This transition 
would be of the nature of a phase transition, like 
condensation, for the gas. Accordingly we could 
say that a transition will take place when the 
total energy is zero. This last hypothesis com. 
bined with the usual statistical method leads to 
results which can only partially be checked by 
experiment. One has to make additional as 
sumptions in order to obtain explicit expressions 
for the transition temperatures. These computa 
tions have been made for a number of gases, and 
a good agreement with experiment is found. 


2. ANALYSIS OF M(c) 


We consider a gas of N particles (x;, yi, 2), 
(¢=1,2,---N) in xyz space. All masses are 
assumed to be equal to m. The distance between 
particle < and particle 7 is denoted by r;;. With 
the mutual interaction between the molecule 
given by f(r)=ar-*—br-™ the total potential 
energy is 
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MECHANICS OF THE CONDENSATION OF GASES 


The energy integral is given by 


(2) 


*Zy) =C. 


2 = (p22 + pui?t+ pz”) + V(x, °° 
2m ; 


Here pzi=m/(dx;/dt), etc. Equation (2) defines the 
locus M(c) in the 6N-dimensional space of the 
variables (x1, ---2v, Pri, +++ pen). 

Since the p’s appear here only as the sum of 
their squares, it follows that the structure of the 
locus M(c) is determined by the structure of the 
locus Vc in the space of (x1, - - -2v). The word 
structure is meant here in the sense of topology, 
i.e., of properties invariant under a continuous 
one-to-one transformation with continuous in- 
verse.-? However, we shall also remark on 
metrical properties. 

To determine the structure of the locus Vc 
we first consider that of V=c, from which the 
former can be deduced. To that end we consider a 
ray from the origin in the (x, - - -2v) space: 


x:=A4it, ¥ix=uil, 2,;= vit, 


D> (AP2+ui?+y,?) =1. 


(3) 


This ray meets the hypersurface V=c in points 
determined by the equation 


(a/t") — (8/t") =c, (4) 


where 


a= > a/pis>0, B=L b/pij>0, 
pig=L(Ai As)? + (ui — we)? + (Hi — vp)? ]. 


It follows that for c=0 there is in general exactly 
one intersection point, while for c<0 there may 
be two points or one or none. There are certain 
exceptional rays which will never meet the 
hypersurface V=c, namely, those for which at 
least one pij=0. These rays correspond to colli- 
sion configurations, configurations in which two 
or more molecules coincide. Because of the energy 
integral, such configurations can never arise for 
any finite value of c. The family of exceptional 
rays forms a point set of 3N—3 dimensions and 
will hence have only a slight effect on the struc- 
ture of the 3N-dimensional locus Vc. 


‘For references to topology, see P. mpeg and H. 
Hopf, Topologie I (Berlin, 1935), especially 1-22. 
* See W. Kaplan, Am. Math. Month. 59, 316.323 (1942). 


3. STRUCTURE OF M(c) FOR c=0 


In this case the above analysis shows that the 
locus VSc has the structure of the exterior of a 
(3N—1)-dimensional sphere (minus certain ex- 
ceptional rays, corresponding again to collision). 
It follows that the energy integral defines a locus 
M(c) which is unbounded in all directions with re- 
spect to the configuration coordinates (x, - - -Zy). 
This type of manifold may be realized as follows: 

Start with the (6N —1)-dimensional sphere 


LD (p27? +puPt+pe*) +L (x2+y92+27)=1. (6) 


Now perform an inversion by reciprocal radii in 
the hypercylinder 


LD (x?+92+2,*) =1. 


That is, we replace each point (x, -- 
-+pay) by the point (x1', ---2y’, pzy’, 
where 


‘ZN, pri, 


ne - pen’), 


pri=per', +++ pav=pay’; 


x1 Xe 2N 
xy’ Xe’ Zy’ 


’ 


X1x1;'=0, - + -tyey’=0; 


(xs? + + + 22) (sey? + + + 8y’2) = 1. 


Finally remove the points on the exceptional rays. 
The result is a locus having the same topological 
structure as M(c) and having the same position® 
in the phase space as M(c). 


4. STRUCTURE OF M(c) FOR c=0 


Here we need a special lemma whose proof is 
omitted : 

Lemma. There is a constant — 6, 6>0, such that 
for —8<c<0 Eg. (4) has exactly two distinct roots 
on every non-exceptional ray. The least possible 
value of —6 is 


—— (:-*)(" il (s 3S aa 


Suppose now —5<c<0. Then VSc holds be- 
tween two hypersurfaces which are topologically 
equivalent to two concentric spheres. (Again ex- 
ceptional rays must be removed.) Thus in this 
case the energy integral defines a locus which is 
essentially bounded. It becomes unbounded only 


* See reference I, p. 2. 
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in those directions which correspond to the ex- 
ceptional rays. A further analysis shows that the 
structure of the manifold is the same as that of 
the topological product of a 3N-dimensional 
sphere and a (3 —1)-dimensional sphere (minus 
certain exceptional points); that is, the same as 
the locus 

3N+1 


3N 

>> u;*=1, > v?= 

i=1 i=1 
in the space of the variables (u, ---ugy, 1, 
++ -Ysv41). For cS=—6, M(c) will collapse even 
further and shrink within itself until for c= —e, 
where — eisacertain constant not yet determined, 
—e>-—, M(c) disappears. 

The value —4 is a second critical value for the 
manifold M(c), but its physical significance has 
not been completely determined. Mayer and 
Mayer‘ in their statistical theory of condensation 
also obtain a second characteristic temperature. 
It is possible that our second change for c= — 64 is 
related to their T,,. 


5. c=0 AS A CRITICAL VALUE FOR THE ENERGY 


The fundamental hypothesis introduced is that 
when the total energy of a gas changes from a 
small positive value to a small negative value a 
transition can occur. It is further assumed that 
this transition is the condensation of the gas. On 
the basis of the model assumed and these as- 
sumptions we can give the following description 
of the mechanism of condensation. 

The molecules are pictured as mass points. 
Because of the nature of the potential energy 
function, a collision in the strict sense, i.e., some 
r;3=0, is impossible. (Collisions in the sense of 
molecules coming so close together that they 
exert considerable forces upon each other will of 
course occur. In statistical mechanics one some- 
times represents such a collision by a discontinu- 
ous collision between non-interacting elastic 
spheres, of finite radius; this is useful for compu- 
tations, but does not represent the actual physical 
state of a collision.) This means that, if the walls 
of the containing vessel were not present, the 
motion of the gas particles would be completely 


4J. E. Mayer and M. Goeppert-Mayer, Statistical Me- 
chanics (John Wiley & Sons, Inc., New York, 1940), p. 308. 
lig and S. F. Harrison, J. Chem. Phys. 6, 87 


without singularity; any motion would be de. 
termined by analytic functions for all time. The 
motion of the system would be represented by a 
smooth curve on the energy surface and it would 
be natural to correlate the properties of the 
system with the properties of this curve as a 
whole. 

In order to take the walls into account, we 
notice first that collisions with the walls do not 
change the energy of the colliding molecules, and 
hence the total energy remains unchanged. The 
only effect of the walls is a discontinuous one: 
Every particle on reaching the wall changes its 
direction continuously but conserves its energy, 
Since the collisions with the walls are isolated 
events, between collisions the system will move 
with regularity as before. On the energy surface 
this means the following: The surface is cut off 
by acertain boundary (corresponding to a certain 
volume) so that the allowed portion is finite in 
extent. (Note that the velocity coordinates are 
always bounded to start with, owing to the energy 
equation.) A trajectory now appears as a curve 
which is smooth until a collision with the wall 
takes place. Thus the trajectory jumps from one 
point of the boundary of M(c) to another point 
and starts from there another continuous path 
until again a collision with the walls takes place, 
and this process is repeated ad infinitum. Thus 
one trajectory appears as a multitude of pieces of 
different trajectories, no two pieces intersecting, 
The effect of the walls may accordingly be 
described as one of continuously varying the 
initial conditions. 

Now the question is how the structure of the 
whole manifold M(c) is related to the character of 
this multitude of trajectories and, in turn, how 
these are related to the structure of the gas. The 
explanation offered is as follows: If c20 the 
surface is infinite in extent in a very strong sense, 
while when c <0 the infiniteness occurs only ina 
small set of directions. It follows, therefore, that 
when c=0 the probability of any orbit’s leaving 
the finite part of the energy surface in a given 
time interval is large. Hence for most initial con- 
ditions the probability of meeting the walls in 
a given time interval is large. If c<0, the orbits 
will wander around for a long interval before 
reaching the boundary. In fact there is a very 
large probability that in a given time interval the 
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orbits remain completely in a bounded part of 
space, and never meet the boundary at all. Thus 
the discontinuous orbit is of a different nature, 
the system shows much less tendency to spread 
out, the gas condenses.® 

It may be remarked that if the walls of the 
vessel are introduced as a force field, instead of 
in the qualitative way as above, the essential 
features of the argument are preserved. The case 
c=0 still appears as a critical value for the 


energy. 
6. APPLICATION OF THE CRITERION 


In order to test the validity of the assumptions 
made, a number of calculations have been made 
to compare the values obtained, on the basis of 
this theory, with the experimental values. In 
order to carry through the computation, one has 
to make some very strong additional assump- 
tions, but it is reasonable to assume that at least 
the order of magnitude should be predicted cor- 
rectly, if the proposed criterion has any value 
whatsoever. 

We assume a gas of N= M? molecules in a cube 
of volume A* and we assume that the molecules 
form, as a first approximation (to be generalized 
below), a cubic lattice placed at the points 
(ha, ka, lx), where h, k, 1=0, +1, --- and 
a=A/M is the nearest distance between mole- 
cules. We use the potential function V as given 
above, but now use numerical values for the 
constants: »=9, m=6, a, 6 as given for various 
gases in Table I. 

The total potential energy of this lattice can 
then be written as 








NfaB bC 
v=5(—-= (8) 
2\a*® a’ 
with 
1 
= 9 C= . 9 
pe (h2+k?+12)%2 p> (h?+k2+12)3 (9) 


Assume now that we have one mole of gas of 
volume v, so that N= No, Avogadro’s number. 
Then, since a= A/M, the potential energy can be 


5 A complete proof of the asserted properties of the orbits 
would require a delicate use of an ergodic-like hypothesis 
or one of metrical transitivity. See G. D. Birkhoff, Proc. 
Nat. Acad. Sci. 17, 650-660 (1931), and J. C. Oxtoby and 
S. M. Ulam, Ann. Math. 42, 874-920 (1941). 
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written as 





N,* /aBN, 
=— -0c). (10) 
2v? v 


The total kinetic energy we assume to be 
3NokT/2; this assumption defines the tempera- 
ture T for this theory. The condition: kinetic 
energy + potential energy = 0 then becomes 











N,? aBN> 
$NokT =—{ 6C— (11) 
2v? v 
or 
1\? 1\3 
T=00(-) -1a(-) (12) 
v v 
where 
NoeC BN, (13) 
——s 


If we determine the numerical values of the con- 
stants o and ¢ assuming a cubic lattice, we find 


o=7.5117X10%, ¢=3.5939108. (14) 


This leads to predictions of the critical tempera- 
ture J, from the known values of the critical 
volume v,. The results are listed in column three 
of Table I. 

Now the assumption of a cubic lattice is not 
justified by experimental evidence. A better as- 
sumption would be that for all gases the spatial 
distribution is a unique one up to a similarity 
transformation. This is really equivalent to as- 
suming the law of corresponding states, which, in 
a sense, is already implicit in this computation 
since we work with a two-parameter potential 
function. It would follow that a formula like (12) 
would still hold, with the same o and ¢ for all 
gases; however, their values are no longer given 
by (14). [This follows by the same reasoning that 
led to Eq. (4). What the exact spatial struc- 
ture is remains to be determined ; the values of 
and ¢ of course depend on that structure.*® It is 
interesting to note that Lennard-Jones and 
Devonshire’ in a series of papers on the theory of 
dense gases and liquids have worked with a 


* In order to compute ¢ and ¢ one must have some infor- 
mation about the distribution of the molecules in a dense 
gas. The actual computations would thus have to be along 
the lines indicated by Mayer and his co-workers. Cf. J. E. 
Mayer and E. Montroll, J. Chem. Phys. 9, 2 (1941). 

7J. E. Lennard-Jones and A. F, Devonshire, Proc. 
Roy. Soc. A163, 53-70 (1937) ; 165, 1-11 (1938) ; 169, 317- 
338 (1939) ; 170, 464-484 (1939). 








20 W. KAPLAN AND M. DRESDEN 


TaBLe I. Compilation of results. Column 3 is based on 
a cubic lattice assumption. Columns 5 and 7 are based on 
the maximum point assumption. 











1 2 3 4 5 6 7 

Te Te Te Ve Ve 
a-10% §«6-108 §~= pred. obs. pred. obs. pred. 
Hydrogen 5.1 17.8 25 33.3 39 64.9 47.6 
Helium a 2.37 44 5.3 9.3 57.8 35.8 
Neon 3.5 14.5 43.6 44.5 45.1 41.7 40.0 
Nitrogen 132 227 145 126 122 90.1 97.0 
Argon 76.8 170 161 151 151 75.2 75.2 
Krypton 193.6 346 219 211 200 107 93.5 

Xenon 583 877 367 290 360 114 111 








similar model and found it to give a consistent 
theory and good agreement with experiments. 
In order to determine the best values of ¢ and 
ft, we first note that according to our theory 
Eq. (12) should give the relationship between 
density and temperature of saturated vapor. 
According to (12) that temperature rises from 
zero to a maximum as the density increases, and 
then decreases. The first part of the curve is thus 
consistent with experiment; the second part has 
not as yet been explained. If the first part of the 
curve is correct, the maximum temperature 
should be 7, and the corresponding volume 2,. 


Using this we find 
4 a} 3fa 
ee ae, v=. (15) 
27 ga? 2ab 


A simple way of determining o and £ would be 
to use the known values of 7, and v, for one 
particular gas to compute o and ¢ by means of 


(15), and then use (15) further to compare the 
critical volumes and temperatures for other gases. 
The results are shown in columns 5 and 7 in 
Table I, with argon as the standard gas. The 
agreement so obtained is in general very satis- 
factory, especially when one considers the un- 
certainty in the force constants. 

It should be remarked, however, that a very 
important factor contributing to this agreement 
is undoubtedly the fact that we use the law of 
corresponding states on the one hand, while we 
use some experimental data from the critical 
region to find o and ¢ on the other.® If we use the 
formula (15) for JT, and v, directly, using for o 
and ¢ the cubic lattice values and for a and b 
values as given from experiment, we have a much 
more a priori theory. If this is done, values are 
obtained having the right order of magnitude; 
they are, however, not so close to the observed 
values as the ones computed in the other way. 
This is to be expected, since in this computation 
the cubic lattice assumptions enter explicitly. 
We can conclude that, in general, the agreement 
with experiment is satisfactory. 

In a paper to follow a more detailed account 
will be given of the method. The connections 
with other theories of condensation will be given, 
and a more refined comparison with experiment 
will be undertaken. 


8 This point was first remarked to us by Professor G. E. 
Uhlenbeck. 
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It was found that reversing the current through the nickel alters the change of magnetism 
both when it is abruptly established and when it is abruptly stopped. In the sixteen different 
cycles considered, invariably a current flowing with the flux opposed magnetization, while a 
current flowing against the flux favored magnetization. A thorough analysis of the numerical 
data shows that in addition to the well-known “shock effect’’ equivalent to a mechanical jar, 
there is another effect which is reversed by reversing the direction of the current with reference 
to the flux. This effect apparently is not due to the circular flux set up by the current nor to 
any dissymmetry in the bar, but instead appears to be some kind of action between the current 
itself and the flux. This hypothesis was strengthened by a supplementary experiment showing 
that the resistance of a nickel bar is decreased when it is strongly magnetized by an alter- 


nating field. 





INTRODUCTION 


HE apparatus used in this investigation was 

identical with that used by Dr. Doolittle 
and myself in our experiments with residual 
magnetism in iron.’ In our published article the 
various circuits and coils have been fully de- 
scribed! and only a brief recapitulation is given 
here. The sample to be tested was a short 
cylinder of nickel, kindly donated by the Inter- 
national Nickel Company. It was cored and 
slotted, so that its section is like a thick and 
nearly closed letter C. The purpose of the slot 
was to eliminate as far as practicable the circular 
flux set up by the longitudinal current through 
the bar. The sample was placed on the axis and 
at the center of a long solenoid in which a current 
of about twoamperes sets up a strong magnetizing 
field. Surrounding the solenoid is a coil of many 
turns of fine wire connected to a ballistic galva- 
nometer. The magnetizing current could be 
either made-or broken abruptly by means of a 
switch, or gradually by moving sliders along the 
coils of three drum rheostats in series, thus 
varying the current from 0.0004 to 1.8 amperes in 
a few seconds. 

The longitudinal current, to be here denoted 
simply by J, was always made or broken abruptly, 
and its effect, as Dr. Doolittle and I had previ- 
ously observed, was equivalent to a mechanical 
shock causing a sudden increase or decrease in the 
bar’s magnetism according to circumstances. 


1H. A. Perkins and H. D. Doolittle, Phys. Rev. 60, 811- 
817 (1941). 
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This change in flux caused the galvanometer to 
deflect by an amount dependent on the violence 
of the mechanical blow or on the magnitude of 
the longitudinal current, and of course on the 
magnetic condition of the sample. 


PROCEDURE 


In the present investigation the value of J was 
always two amperes, and the magnetizing current 
was always brought, either quickly or slowly, toa 
final maximum very close to 1.8 amperes. The 
procedure involved carrying the nickel cylinder 
through a quasi-hysteresis cycle, and observing 
the deflections of the galvanometer at the 
vertical ends of the cycle as shown in Fig. 1. As 
these deflections are proportional to the change in 
flux caused by making or breaking J, they alone 
are quoted in what follows which is-only con- 
cerned with relative values. 

Sixteen cycles were examined, each of which 
may be represented qualitatively by one of the 
diagrams in Fig. 1. In (a), the nickel cylinder was 
magnetized following a curve similar to a typical 
B—H curve from a to }, while the current J was 
flowing through it. At b, J was stopped. This 
caused an increase in magnetism bringing the 
flux abruptly to c which may be regarded as the 
normal condition for that field intensity. The 
field was then decreased to d when J was re- 
established and the condition at a was recovered 
by a decrease of flux proportional to da. The 
increase from } to c and the decrease from d to a 
indicate a hysteretic lag which always occurs in 
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nickel, because the property of ‘‘overshooting the 
mark”’ previously found in Norway iron! does not 
exist with nickel or cast iron. 

It should be noted that Figs. 1 are misleading 
in one respect. They exaggerate the deflection of 
the galvanometer at the two ends of the cycle. 
Although these ‘‘throws’’ were often over 100 
scale divisions, the reading that would have been 
produced by the abrupt magnetization or de- 
magnetization of the cylinder would have run to 
a thousand or more. An accurate drawing of such 
a cycle would show two curved lines very close 
together. 

Another point to be noted is the absolute 
necessity ef going through each cycle four or five 
times before any readings are made, in order to 
get the nickel accustomed to the regime. Nickel, 
and iron even more so, have amazing memories, 
and it takes considerable treatment to make them 
forget their immediate past. Taking this pre- 
caution, a given cycle in the case of nickel could 
be repeated again and again, or recovered later 
on with surprising consistency. Unfortunately 
Norway iron did not behave so well, although 
some of the same effects as with nickel were 
observable. 

The sixteen cycles were obtained from the 
various possible combinations of the following 
factors. The magnetizing field may be established 
or discontinued either abruptly or slowly. The 
current may be flowing while the flux is in- 
creasing and not flowing while it is decreasing, or 
flowing during decreasing flux and not flowing 
during increasing flux. It may also flow during 
both processes or not flow during either except at 
the ends of the cycle where a momentary opening 
or closing of the switch is needed to obtain the 
galvanometer throw indicating the change of 


Fic. 1. 
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flux. As there are two such abrupt changes for 
each cycle, indicated by bc and da in Figs. 1, the 
16 cycles give 32 different readings of the galva- 
nometer. Then by reversing J, we obtain 32 more 
making a total of 64 distinct observations. Each 
of these observations was the average of four 
similar readings which never varied more than 
about seven percent from the average and 
generally very much less. 

As a result of a careful study of the 64 different 
cases observed, it was found that, without ex- 
ception, in stopping the longitudinal current J 
after magnetization, or in starting it after 
demagnetization, the galvanometer throw is 
greater when IJ is with the flux (denoted by +J) 
than when it is against it (denoted by —J). In 
starting J after magnetization or stopping it after 
demagnetization, the galvanometer throw is less 
when J is positive, and greater when it is nega- 
tive. The meaning of these differences may be 
understood by referring to Figs. 1. In (a) the 
deflection due to the flux change b’c when —J 
flowed during magnetization is less than that due 
to bc when +I was flowing. Thus —J favored the 
growth of the flux while +J held it back to the 
lower point b. 

Figure 1(b) shows what happens at the other 
end of the cycle when J flowed during the process 
of demagnetization. When the current is with the 
flux (+J) the curve descends to d, while with —J 
it only goes to d’. Thus there is more residual 
magnetism when the current is against the flux 
than when it is with it. So here again —J favors 
the magnetic state as compared to +J which 
presumably opposes it. 

Figures 1(c) and (d) illustrate abruptly starting 
I at the two ends of the cycle. In (c) starting —J 
causes the larger throw (0 to c’) indicating a 
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INTERACTION BETWEEN CURRENT AND FLUX 


TABLE I. Results for eight different cases. 








Reading when J 
flows in nickel 


Procedure +I —I 


Relative decrease 
of deflection 
when J is 
reversed 


Reading when J 
flows in axial wire 


+1 of 


Relative decrease 
of deflection 
when J is 





+64 +50 
+51 +58 
+43 +34 
+31 +40 
—100 —140 
—146 —105 
—-23 -—49 
—51 —25 


. Break after slow magnetization 
Make after slow magnetization 
Break after quick magnetization 
Make after quick magnetization 
Break after slow demagnetization 
Make after slow demagnetization 
Break after quick demagnetization 
. Make after quick demagnetization 


FOP A > Po 


+66 +64 
+59 +64 
+34 +51 
+49 +54 
—113 —120 
—120 —112 
—47 —45 
—48 —48 








greater increase of magnetism than that caused 
by starting +J which carries the flux only to c. 
In (d) starting —J after demagnetization brings 
the flux from d to a’, while starting +J brings it 
to a. But a’ indicates more residual magnetism 
than a, so in this fourth case also, —J favors the 
magnetic state while +J tends to diminish it. 
The average differences expressed in galva- 
nometer scale divisions between the effects of +J 
and —J for the four cases just described and 
based on all 64 observations are as follows: 


1. Break after magnetization: 12 mm due to bb’, 
Fig. 1(a); 

2. Break after demagnetization: 33 mm due to 
dd’, Fig. 1(b); 

3. Make after magnetization: 8 mm due to cc’, 
Fig. 1(c); 

4. Make after demagnetization: 34 mm due to 
aa’, Fig. 1(d). 


It is evident that the phenomenon is much 
more pronounced after’ demagnetization than 
after magnetization. This is to be expected, be- 
cause as we approach saturation the magnetic 
state becomes increasingly rigid. Also the relative 
changes of flux are found to be larger after quick 
than after slow demagnetization, although the 
actual changes involved are smaller, because 
quick demagnetization results in a smaller resid- 
ual flux than when the magnetism is gradually 
reduced. 

In Table I are listed eight essentially different 
cases. Each of the galvanometer throws indicated 
under +J or —IZ is the average obtained from 
four different cycles. These four values agree 
fairly well among themselves, especially after 
demagnetization, which shows that an individual 
case is not much influenced by the rest of the cycle. 


DISCUSSION OF RESULTS 


It will be seen that after magnetization the 
deflection is listed as positive to indicate an in- 
crease of flux, while the negative sign after 
demagnetization indicates a decrease as would be 
expected. The first two columns of figures are 
those already discussed. The other values for an 
axial current will be explained later. The column 
indicating the relative change due to a reversal of 
I shows that in some cases the change is very 
marked and is never less than 12 percent. This 
effect can be explained by supposing two major 
effects, A and B, and a third minor one, C. The 
first of these, A, is the well-known “shock effect”’ 
which alters the magnetization as a result of a 
mechanical jar or a sudden rush of current, inde- 
pendent of its direction. The second, ‘‘B effect,” 
accounts for the change of deflection when J is 
reversed. It is this phenomenon with which the 
present discussion is chiefly concerned, and it is 
apparently not accounted for by classical electro- 
magnetic theory. . 

The C effect is due to the action of the circular 
flux set up in the bar by the current flowing 
through it. This circular flux might be supposed 
to explain the effect of reversing J, but it was 
shown to be much too small to account for the 
magnitude of the observed phenomena. A 
crucial test proving this fact was made as follows. 
Instead of sending J through the nickel, it was 
carried in a wire along the axis of the cored 
cylinder. Thus only two effects were possible, 
namely, shock and the interaction between the 
circular flux set up by the axial current and the 
longitudinal flux in the bar. This interaction, in 
spite of the slot, designed to minimize circular 
flux, must also occur to some extent when the 
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current flows through the metal. But calculation, 
as stated in the earlier article,! shows that the 
circular flux is 2.5 times as large when the current 
is axial, so we should now expect it to assume 
some importance. Thus though the shock effect 
must be independent of the direction of J, the 
influence of the circular flux should result in 
somewhat different galvanometer throws when I 
is reversed. 

The results obtained with the axial current, 
keeping all other conditions the same as before, 
are given in the last three columns of Table I. It 
will be seen that the galvanometer throws are 
similar to those previously obtained. But, though 
of the same general magnitude, the lower value 
of the various pairs is now nearly equal to the 
higher values and the proportional decrease does 
not exceed 9 percent obtained in case number 4 
as compared to 23 percent previously observed. 
It goes down practically to 0 percent as compared 
to 51 percent in case number 8. The differences 
caused by reversing J when it is axial are proba- 
bly due to circular flux, but though observable in 
all but case number 8, it is relatively small, and 
since the axial current is 2.5 times as effective in 
setting up a circular flux as the current through 
the bar itself, we may be justified in dividing the 
percentages by that figure. Then the largest pro- 
portional change of 9 percent would be reduced to 
3.6 percent as compared to 23 percent when the 
B effect as well as the C effect is present. Thus the 
C effect is shown to contribute something but not 
a great deal to the change of deflection brought 
about by reversing J. 


ANALYSIS OF RESULTS 


Since we are now justified in neglecting the 
relatively small C effect as an explanation of the 
observed phenomena, there remain the A and B 
effects to be considered. These may be separated 
as follows. The shock effect in nickel always in- 
volves an increased flux after magnetization and 
a decreased flux after demagnetization, for nickel, 
unlike Norway iron, shows no evidence of “‘over- 
shooting the mark.’’ This means that shock tends 
to carry the flux nearer to its final stable value. 
But the B effect tends to help or oppose this re- 
turn to a stable state according to its direction, 
and whether or not reaching the stable state calls 
for an increase or decrease of magnetization. 


j 





PERKINS 


Thus if we call the change of flux caused by shock 
AA, and that caused by the B effect AB, we have 
AA+AB=M when the effects are in the same 
sense, and AA —AB=N when they are opposed. 
Thus AA=}3(M+N) and AB=}(M—N). Ap- 
plying this method to case number 7 (the most 
conspicuous one) we obtain from the galva- 
nometer readings (proportional to AA and AB) 
AA = }(23+49) = 36 which is proportional to the 
shock effect, and AB=4(49—23)=13 which re- 
presents the B effect caused by +J flowing during 
the process of demagnetization. The same analysis 
applied to case number 3, for instance, gives 
AA = 38.5 and AB=4-.5. If this analysis is applied 
to the corresponding cases when the current is 
axial, we may separate the A and C effects. Thus 
case number 7 gives AA = 46 and AC=1, while in 
case number 3 AA=52.5 and AC=1.5. These 
results seem to justify ignoring the effect of 
circular flux in estimating the B effect when the 
current flows through the bar, although strictly 
speaking the values 13 and 4.5 obtained above 
are partly due to C as well as B but with AC 
considerably smaller than 1 and 1.5 for reasons 
already mentioned. Aside from circular flux, 
other possible explanations of AB caused by re- 
versing J are dissymmetry of the bar itself due to 
crystalline structure, and the effect of the mag- 
netic field of the earth to which it was parallel. 
Both these were shown to be negligible by re- 
versing the magnetizing current. When this is 
done, substantially the same values as before are 
recovered provided J is also reversed. Several 
cycles were carried through in this way, with 
practically the same galvanometer readings as 
before. 


DISCUSSION OF RESULTS 


To sum up the results, it seems safe to conclude 
that the B effect represents some apparently new 
kind of interaction between current and flux ®, 
favoring ® when it flows in the opposite direction 
and opposing ® when in the same direction. Thus 
we may imagine a pure shock effect caused as by 
a mechanical jar tending toward magnetic sta- 
bility with a bias due to B added to it or sub- 
tracted from it according to the direction of the 
current with respect to the flux. This effect exists 
in nickel, but to a much smaller degree, if at all, 
in a bar of Norway iron which was carefully ex- 
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amined with a view to detecting a possible 
alteration in the galvanometer deflection when J 
is reversed. 

Finally the question naturally arises: Is there a 
reciprocal phenomenon? That is, if making or 
breaking a longitudinal current, aside from the 
effect of shock, alters the magnetic state of the 
iron, does a varying magnetic field alter the cur- 
rent? This would mean a change in the effective 
resistance of the bar. To test this possibility, bars 
of soft iron and of very pure nickel kindly 
donated by the International Nickel Company 
were wound with a double-layer solenoid and 
magnetized by an alternating current of 60 cycles. 
Their resistances both magnetized and demag- 
netized were measured in a Kelvin double bridge. 
The iron exhibited no observable change of re- 
sistance, but the resistance of the nickel was 
decreased by more than 0.6 percent when the 
field was applied, and remained at the lower 
value as long as the alternating current was 
flowing. There was of course a decided gradual 
increase of resistance caused by rising tempera- 
ture, but during the short time needed to meas- 


ure its new and smaller value, when the alter- 
nating field was applied, the rise due to higher 
temperature was much too small to obscure the 
phenomenon. 

It has long been known that a steady longi- 
tudinal field raises the resistance of nickel meas- 
ured in the same direction, that the relative 
increase, AR/R, is unaltered by reversing the 
field, and that this magnetoresistance effect is 
associated with the still older phenomenon of 
magnetostriction.? But apparently the effect of a 
rapidly alternating field has not been examined, 
and it is certainly surprising to find that such a 
field lowers the resistance even more than it is 
increased by a steady field. This increase was also 
measured with similar excitation by the author 
and it was found to be between 0.5 and 0.6 
percent instead of over 0.6 percent with an 
alternating current. In both cases the flux density 
was about 1900 gauss which was unnecessarily 
high, as the saturation flux density for this effect 
is of the order of 100 gauss. 


?L. W. McKeehan, Phys. Rev. 36, 948 (1930). 
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E have performed a systematic study of 

the electret effect,! guided by the idea 
that it must be related to the various other 
aspects of anomalous behavior, of solid dielec- 
trics. The manufacture of electrets, as it is usually 
accomplished, involves a heat treatment. Prior, 
therefore, to the investigation of the fields of elec- 
trets, we studied the influence of the temperature 
on the dielectric absorption. A preliminary report 
of our results is presented here. 














TERMINOLOGY 





Following Gemant,' we call homocharges the 
true (surface or space) charges, the sign of which 
is the same as that of the corresponding polariz- 
ing electrode, and we call heterocharges the true 
charges of contrary sign. Both kinds of charges 
coexist in a polarized dielectric. In the shorted 
capacitor, their superposed fields induce on the 
electrode a charge g, which is measured by 
lifting the electrode. An apparent surface 
charge —q is then ascribed to the dielectric. 










INFLUENCE OF TEMPERATURE ON 
DIELECTRIC ABSORPTION 







Here we are not concerned with the measure- 
ment of isothermic current—time curves, but 
with non-isothermic effects, which have a bearing 
on the electret effect. Their existence, suggested 
by former papers,” could be verified : A capacitor 
containing an absorptive dielectric (like carnauba 
wax) is charged at a high temperature and then, 
always connected with the voltage source, cooled 
until a considerably lower temperature is 
reached; finally it is short-circuited. Neverthe- 
less, only a partial discharge occurs. A part of 















1Cf. M. Eguchi, Phil. Mag. 49, 178 (1925). S. Mikola, 
Zeits. f. Physik 32, 476 (1925). E. P. Adams, J. Frank. 
Inst. 204, 469 (1927). M. Ewing, Phys. Rev. 36, 378 (1930). 
R. D. Bennett, Phys. Rev. 37, 103 (1931). A. Gemant, 
Phil. Mag. 20, 929 (1935). G. Groetzinger and H. Kretsch, 
Zeits. f. PI ysik 103, 337 (1936). P. A. Thiessen, A. Winkel, 
and K. Herrmann, Physik. Zeits. 37, 511 (1936). W. M. 
Good and J. D. Stranathan, Phys. Rev. 56, 810 (1939). 
B. Gross, Anais Acad. Bras. Ci. 15, 63 (1943). 

2H. Frei and G. Groetzinger, Physik. Zeits. 37, 720 
(1936). O. G. v. Altheim, Ann. d. Physik 35, 417 (1939). 
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the quantity of electricity stored up during the 
preceding charging cycle behaves like a ‘‘frozen”’ 
charge and is retained within the dielectric as 
long as the temperature is kept low. The value 
of the frozen charge depends on the nature of 
the dielectric, on the charging time, and on the 
difference of temperature between charge and 
discharge, but it is independent of the rate of 
cooling. The explanation of the effect is found 
in the enormous increase, with rising tem- 
perature, of both the charging and the dis- 
charging velocity. The true absorption capaci- 
tance of capacitors with solid dielectrics is fre- 
quently much higher than is usually believed. 
But the charging velocity at room temperature 
is so low that at such temperatures the capacitors 
are never charged completely. On the contrary, 
in the first stage of our experiment the capacitor 
is almost completely charged in a brief time, 
because it is charged at a high speed at the high 
temperature. In the second stage (beginning with 
the short circuit), the capacitor is far from being 
completely discharged even after a long time, 
because at the low temperature it is discharged 
at a low speed; thus the dissipation of the ab- 
sorbed quantity of electricity is so small that a 
part of it is trapped in the dielectric. Whatever 
may be the nature of these frozen charges, they 
must produce electric fields. This type of dielec- 
tric absorption yields therefore necessarily the 
electret effect ; but it does not explain everything, 
because (a) the mechanism, which according to 
current theories (dipole and ionic theories) is 
responsible for dielectric absorption, accounts 
only for the heterocharges, not for the homo- 
charges of electrets; and (b) electrets can be 
made without heat treatment. 


DECOMPOSITION OF THE CURRENT 
INTO DISPLACEMENT CURRENT 
AND CONDUCTION CURRENT 


We have investigated the nature of the 
homocharges by varying the intensity and time 
of application of the polarizing field, by modi- 
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EXPERIMENTS 


fying the temperature, and by speeding up the 
disappearance of the field of the electret by 
heating. During the discharge, the total current 
J and the charge g of one of the electrodes were 
measured. According to the equation J =i+dq/dt, 
this method permits us to split up the current 
into its two components, the conduction current 
4 and the displacement current dg/dt, both 
referred to the interface dielectric-electrode. The 
conduction current originated by the decay of a 
homocharge has the usual direction of the dis- 
charging current, but the corresponding displace- 
ment current flows in the direction of the charg- 
ing current; just opposite in direction are the 
currents for heterocharges. The decay of the 
homocharges proceeds partially by conduction 
across the interface dielectric-electrodes, par- 
tially by conduction within the dielectric. If the 
second effect is sufficiently strong, the total dis- 
charging current may reverse.* These current 
reversals are not to be confused with those which 
are accounted for by the principle of superposition 
and which are observed after a brief application 
of a voltage of one sign succeeding a prolonged 
application of a voltage of the same value but 
contrary sign.‘ 


MEASUREMENTS ON CURRENTS AND 
CHARGES IN ELECTRETS 


Even without heat treatment, numerous die- 
lectrics produce strong fields after the application 
of high voltages (cf. S. Mikola).1 An important 
factor is the duration of the polarization. The 
results are particularly clear for carnauba wax. 
(1) If the time of polarization has been very short, 
the apparent charge—with the sign of a homo- 
charge—decays at first rapidly. The decay is 
associated with a strong displacement current 
so that the total current flows in the same direc- 
tion as during the charging period. Later the 
apparent charge becomes more stable and the 
current goes down. When the temperature is 
increased, the charge disappears completely and 
a transient displacement current is raised. (2) If 
the time of polarization has been long, shortly 


* Reversals of a similar kind seem to have been observed 
by H. H. Race, Trans. A. I. E. E. 47 (October, 1928) and 
by Sadakichi Shimizu, Tohoku Univ., Sci. Reports 21, 29 
(1932), see especially p. 48. 

*Cf. V. Karapetoff, Trans. A. I. E. E. 45, 124 (1926). 
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after the short circuit an increase, and later a 
maximum of the apparent charge (which again 
has the sign of a homocharge), are observed. The 
current is now mainly conduction current. It is 
the usual discharging current familiar from the 
numerous experiments on anomalous currents in 
solid dielectrics. A later increase of the tem- 
perature is followed by an immediate decay of 
the apparent charge, accompanied, as in the 
first case, by a displacement current. 

The usual heat treatment produces the same 
effect as an enormous prolongment of the 
duration of the polarization. The value of the 
temperature itself has only quantitative in- 
fluence. It need not be raised above the melting 
point of the substance. Strong permanent elec- 
trets can be obtained by polarization in the solid 
state. It becomes now evident that immediately 
after the short circuit the sign of the apparent 
charge is always that of a heterocharge; but 
sooner or later a change of polarity occurs and 
it turns into that of a homocharge. Simul- 
taneously a strong discharging current, almost 
entirely conduction current, is observed. If the 
temperature of the shorted system is raised 
temporarily, the value of the apparent charge 
(homocharge) is greatly enhanced, but it returns 
slowly to its former level after room temperature 
has been reached again. If the temperature is 
raised and maintained high for a long time, the 
apparent charge disappears after a transitory 
increase and a strong discharging current surges, 
again a conduction current. 

The systematic deviations from the principle 
of superposition, which at high voltages become 
obvious, are accounted for by the presence of 
homocharges. 

Ionization in an air gap eventually existing 
between dielectric and electrode is not sig- 
nificant. We have observed current reversals 
indicating the existence of homocharges on hard 
rubber disks with silvered surfaces. 


THEORY OF ELECTRET BEHAVIOR 


Two mechanisms act simultaneously: (1) die- 
lectric absorption related to the movement of 
ions or the orientation of dipoles in the interior 
of the dielectric originates heterocharges ; (2) con- 
duction in the interphase dielectric-electrode 
which originates homocharges. Without heat 
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treatment, the heterocharges, with a smaller 
relaxation time, disappear spontaneously after 
the system has been shorted, and a pure homo- 
harge remains. The velocity of the transition de- 
pends on the degree of dielectric absorption, 
which in its turn is determined by the duration 
of the polarization. After a heat treatment, a 
considerable fraction of the heterocharge is pre- 
served besides the homocharge in the shorted 
system. The resulting field depends on the values 
and spacial distribution of both charges; its tran- 
sitory increase during a later heating is another 
consequence of the smaller relaxation time of 
the heterocharges. 

In polar substances, the formation of hetero- 
charges is due mainly to dipole orientation. 
Proportionally to the value of these charges 
the electric field increases in the interphase 
dielectric-electrode. With sufficiently high field 
strengths, conduction currents surge in the 
> interphase and, consequently, ions (or electrons) 

are fed into the dielectric or extracted from its 
surface and transferred to the electrodes. Ac- 
cordingly, homocharges appear, initially in the 
form of surface charges, later spraying over a 
certain depth within the dielectric.' This process, 
resulting in a weakening of the field, reduces the 
conduction currents. With the short circuit, a 
part of the dielectric polarization, in phase with 
the electric field, disappears immediately; the 






5This picture is very similar to that given by M. 
Faraday, Experimental Researches in Electricity, Volume I 
(Richard and John Edward Taylor, London, 1839), who 
has already known the existence of homoch : “The 
effects appear to be due to an actual penetration of the 
charge to some distance within the electric, at each of its two 
surfaces, by what we call conduction; so that, to use the 
ordinary phrase, the electric forces, sustaining the induc- 
tion, are not upon the metallic surfaces only, but upon and 
within the dielectric also, extending to a smaller or greater 
depth from the metal linings.” (Art. 1245, p. 390.) 
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rest follows more or less slowly. In consequence 
of this decay of the heterocharge, the homocharge 
component begins to prevail in the resulting 
field. Again conduction currents play a role, now 
carrying away part of the homocharge. It must 
be kept in mind that there is no recombination 
between the two types of charges. Thus, the 
homocharge is released in the same measure as 
the heterocharge decays. If a part of the hetero- 
charge is frozen in, the relaxation time for its 
decay is so long that the homocharge in its turn 
is conserved over extremely long periods of time. 
Such an electret possesses really a volume 
polarization like a magnetized iron bar, but the 
field becomes disturbed by conduction effects 
non-existing in magnetism. If an electret with 
strong homocharges is cut into two pieces of 
about equal dimensions, both surfaces of the 
same piece should show charges of the same sign. 
But conduction currents surge at the newly 
created surface in contact with metal, the in- 
ternal field of the electret now playing the same 
role as previously did the external field. So a 
change of polarity takes place and each piece 
becomes a dipole. 

The effect, which manifests itself in the ap- 
pearance of homocharges, must interfere in 
many problems of dielectric behavior, like 
increase of dielectric loss at high voltages and 
rupture. Its closer investigation may reveal 
further interesting features. 
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On a Recent Paper by Jaffé 


J. R. Prerce 
Bell Telephone Laboratories, Inc., New York, New York 
June 22, 1944 


HAVE been interested in comparing a paper, ‘On the 

currents carried by electrons of uniform velocity,” by 
George Jaffé,! with an earlier paper on the subject by C. E. 
Fay, A. L. Samuel, and W. Shockley.? These treatments 
agree in all cases in which the current is not “space-charge 
limited.” The difference may be illustrated by the following 
special case. 

Suppose we have a thermionic cathode, a positive grid 
close to it, and an anode at the same potential as the grid, 
spaced many times as far from the grid as the grid is from 
the cathode. Imagine that we gradually increase the tem- 
perature of the cathode, and hence the current from it, 
until the cathode current finally becomes space-charge 
limited. 

According to Fay, Samuel, and Shockley, the entire 
cathode current can be transmitted through the grid-anode 
region until a “limiting current” is reached. For cathode 
currents greater than this, and also as an alternative possi- 
bility for currents somewhat smaller, they would regard as 
valid an expression which involves a potential minimum at 
cathode potential between grid and anode, with part of the 
electron flow turned back toward the cathode at this po- 
tential minimum. Under this condition the anode current is 
considerably less than the “limiting current.” Clearly, Fay, 
Samuel, and Shockley regard this condition only as repre- 
senting approximately the actual solution with thermal 
velocities, for they say, “The assumption of equal energy 
electrons leads to indeterminacy at planes where the 
potential and its gradient are both zero. At these planes the 
existence of non-uniform electron velocities will be recog- 
nized insofar as it provides a selective mechanism to resolve 
the mathematical indeterminacy.” 

On the other hand, Jaffé rejects completely any expres- 
sions implying zero electron velocity and infinite charge 
density at a plane, and according to his treatment the final 
value of the anode current is the “limiting current” of Fay, 
Samuel, and Shockley, mentioned above. 

Jaffé says, ‘Now, the velocity of the electrons naturally 
becomes zero at a point where 7 = —1 [see (3) ], and there- 
fore the density m would have to be infinite. We are of the 
opinion that no solution should be admitted which makes a 
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physically measurable quantity infinite. We, therefore, 
conclude that the limiting case X; =0 (on branch IT) should 
be excluded even as a possible form of unstable discharge.” 
This, of course, excludes all solutions for which a part of the 
current is turned back, and indeed, Jaffé says in his intro- 
duction, ‘In general, the minimum will not be low enough 
to stop electrons which have entered the discharge space. 
In spite of this fact the current will be space-charge limited 
(for sufficiently low values of the potential) because the 
solution of the problem does not exist for any given value 
of the plate distance unless the current is below a critical 
value. Thus the space-charge limitation is caused, in this 
case, not by the fact that some electrons return to the 
cathode, but by the fact that only a limited number of 
them can get into the discharge space.” 

It is pretty clear experimentally that electrons can be 
turned back in a space between approximately plane 
electrodes. The writer believes that expressions given by 
Fay, Samuel, and Shockley which involve a plane of zero 
velocity and infinite charge density (the charge is not 
infinite) are physically acceptable as representing a limit 
which the correct solutions involving thermal velocities 
must approach as the parameter expressing the ratio of the 
thermal energy to the energy at the entrance or exit plane 
is made very small. 

In regard to stability, Jaffé’s argument that in cases in 
which the velocity does not go to zero, the solution with 
the lower potential minimum (Fay, Samuel, and Shockley 
C overlap solution) is unstable because “it would require 
work to change the discharge with the higher potential 
curve to that with the lower potential curve” is not deci- 
sive, as by this criterion a bottle is unstable when standing 
upright. Fay, Samuel, and Shockley also regard the C 
overlap solution as unstable, and there are reasons to 
suspect that it is completely unstable. 

In a paper, ““Vacuum-tube networks,”’ by F. B. Llewellyn 
and L. C. Peterson,’ the a.c. problem is treated. For an 
undisturbed current injected into a diode, the relation 
between a.c. (or perturbation) voltage V across the diode 
and a.c. (or perturbation) current J in the diode is given by 
an expression which for very low frequencies reduces to 


V=1(ue+uo)(T/e)(1—$+$8/3)/8. 


Here uw. and wm are velocities at the entrance and exit 
planes, « is the dielectric constant of vacuum, T is transit 
time across the diode, ¢ is a “‘space-charge factor” which for 
the C overlap solution is always greater than unity, and 
8=iwT. We note that for certain small values of 8 such that 


current oscillations can occur in a short-circuited diode 
(V=0). We also note that for ¢ a little less than unity the 
oscillation starts as an exponential decay of perturbation 
current with time, while for ¢ a little larger than unity the 
oscillation starts as an exponential increase of perturbation 
current with time, indicating an instability. 

1 George Jaffé, Phys. Rev. 65, 91-98 (1944). 

2C. E. Fay, A. L. Samuel, and W. Shockley, Bell Sys. Tech. J. 17, 49- 
79 (1938). The same matter was treated by B. Salzberg and A. V. Haeff 


in the ar Review (Jan. 1938). 
3 F, B. Llewellyn and L. C. Peterson, Proc. I.R.E. 32, 144-166 (1944). 
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On the Currents Carried by Electrons of Uniform 
Initial Velocity 
GEORGE JAFFE 


Louisiana State University, Baton Rouge, Louisiana 
July 3, 1944 


N a Letter to the Editor, Mr. J. R. Pierce has compared 
my recent paper under the above title' with an earlier 
treatment of the same subject by C. E. Fay, A. L. Samuel, 
and W. Shockley.? Though Mr. Pierce’s representation of 
my point of view is entirely correct, I should like to add a 
few remarks since it seems to me that a point of principle is 
involved. 

The paper of the three authors was not known to me until 
Mr. Pierce drew my attention to it. However, the treat- 
ment therein is based on the treatment of Langmuir and 
Compton? to which I have referred. 

The original derivation of the Child-Langmuir formula 
was subject to the following two criticisms: 

(a) It does not contain the slightest indication of the 
phenomenon of saturation. 

(b) It makes the solution singular in a mathematical 
sense by admitting an infinite value of the space charge 
density (at the cathode). 

Admittedly the treatment of L. and C. and of the later 
authors, in the case of constant initial velocity, is open to 
objection (b) as regards the plane where the potential 
minimum occurs. It was the only aim in my investigation 
to do away with the mentioned objections. By introducing 
a modified form of boundary condition, I obtained solutions 
which admit of no indeterminacy, which are regular every- 
where (also in the case of space-charge limited currents), 
and which automatically change into the saturated solution 
when a critical potential is reached. 

If Mr. Pierce and others are willing to admit infinite 
values of the space-charge density, there is nothing to be 
said against the treatment of L. and C. except that a very 
puzzling ambiguity arises. Which solution is to be accepted 
as the valid one, the regular or the singular one? From 
general physical and mathematical principles there can 
hardly be a doubt as to the answer. However, there is a 
further argument in favor of the regular solution. 

When the space-charge limitation sets in, the minimum 
of the potential is not yet as low as its lowest possible value. 
(Thus in the case §)= 10 represented by Fig. 3 of my paper, 
the limitation sets in for 7min= —0.45, and for real cases of 
thermionic emission the minimum would be still less low.) 
If now the limitation sets in (e.g., by an increase in the 
potential) according to the solution of L. and C., the mini- 
mum would have to jump suddenly to nmin= —1. It seems 
to me that such a behavior is quite inadmissible and that, 
therefore, the regular solution has to be continued into the 
domain of space-charge limitation. This is exactly what my 
modified boundary condition achieves. 

How nearly the singular solution approximates the 
regular one remains to be discussed. It is easy to show that 
the treatment of L. and C. approximates the regular solu- 
tion in two limiting cases, namely, either (in the notation 
of my paper) when E,/E-~0, or when £o—>5/3, which means 
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when 16mi/(voEo)4x-—>5/3. In the former case the minimum 
will be slight (mmin—>0) and near the cathode; in the latter 
case it will be near the lowest possible value (nmia—> — 1) and 
near the anode. 

In all other cases deviations must be anticipated which 
concern the potential distribution as well as the charac- 
teristic. It will be sufficient to show the discrepancy in the 
latter case. Particularly large deviations are likely to occur 
for relatively large initial velocities, i.e., in the limiting case 
E,/E>>1. In this limit my formula (30) yields the value 
j=8£,!. On the other hand, it can be deduced from L. and 
C.’s formulae (293) to (295) that 7 = Eo!((1+-7)/y)*, where 
v= (22,/i—1)4. Since the saturation value i, is always 
larger than i, y will be a number larger than unity, and the 
factor multiplying Zo! will be between } and } of the value 
in my formula. 

I should like to stress once more that my treatment re- 
gards exclusively the idealization that all electrons have the 
same velocity normal to the cathode. It is only in this case 
that the electrons cannot build up a minimum sufficiently 
low to make them revert (unless a singularity is admitted), 
If there is a velocity distribution, it has been shown long 
ago, and I have stressed the fact in my paper that the 
slower electrons will revert in the field created by the 
swifter ones. 

As to the question of stability, I agree with Mr. Pierce 
that it cannot be answered completely without a special 
investigation. Therefore, I should have called the solution 
with the lower minimum not the “unstable solution” but 
the “‘less stable solution.”” As a matter of fact, I did admit 
a certain stability for the “unstable solution’’—otherwise I 
would not have drawn diagrams giving the potential 
distribution in such cases, For the main point under dis- 
cussion this is quite irrelevant. I pointed out that the 
solutions which make the minimum as low as nmin = — 1 are 
on the unstable branch, but I excluded them as physically 
inadmissible because they make the space-charge density 
infinite. 

1 George Jaffé, Phys. Rev. 65, 91 (1944). 
2C. E. Fay, A. L. Samuel, and W. Shockley, Bell Sys. Tech. J. 17, 49 


(1938). 
31, Langmuir and K. T. Compton, Rev. Mod. Phys. 3, 191 (1931), see 


p. 239, 





Retraction of Stressed Rubber 


B. ADALBERT Mrowca, C.S.C., S. LEONARD Dart, AND EuGENE GuT# 
Department of Physics, University of Notre Dame, Notre Dame, Indiana 
July 3, 1944 


HE process of the retraction of stressed rubber has 
been studied by us for some time, using several ex- 
perimental methods. The most promising method employed 
was the use of a smoked revolving drum driven by a 
synchronous motor. A light stylus was fastened to the free 
end (or any position behind the tip) of a rubber sample; the 
other end of the sample was clamped. The elongation-time 
curves were recorded accurately and directly on the drum 
by the stylus. From these curves the velocity of the tip, as 
a function of elongation (or time), may be obtained by 
differentiation. Figure 1 shows elongation-time curves. 
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Fic. 1. Elongation-time curve for Hevea gum. 


From the curve for about 450 percent initial elongation a 
tip velocity higher than 220 miles/hr. was derived. The 
velocity of the tip is sensibly constant, except for the initial 
acceleration at the beginning of the retraction. This initial 
increase is due to the weight of the stylus and of the rubber 
in the grip (before the sample is released). 

Figure 2 gives the result of a two-stylus experiment, one 
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Fic. 2. Elongation “vs. time with two}styluses, one at the tip and one at 
the middle ‘ofthe samplejofjHevea tread. 


at the tip, the other at the middle. It clearly shows the lag 
of the middle behind the tip and hence the element by 
element nature of the retraction. The middle moves only 
slightly slower than the tip. 

Figure 3 shows the result of a two-stylus experiment for a 
Butyl tread stock. Here, the elongation-time curve for the 
tip has an experimental tail. The lag of the middle behind 
the tip is much smaller. The velocity of the middle point is 
also considerably smaller than that of the tip. The velocity 
of each section of the sample appears to be proportional to 
the distance from the fixed end. 

Figure 4 shows the maximum tip velocities plotted 
against the initial elongation for all these stocks. It is 
noteworthy that GR-S tread, while more sluggish at small 
elongations, compares favorably with Hevea tread above 
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Fic. 3. Elongation vs. time with two styluses, one at the tip and one at 
the middle of the sample of Butyl tread. 


90 percent. The great reinforcement for GR-S and the 
comparatively small one for Butyl are exhibited clearly in 
this graph. 

Hevea gum and Butyl tread are the two extremes. Hevea 
tread, GR-S gum and tread, and Butyl gum were found to 
show an intermediate behavior. 

The present work is the first (to the authors’ knowledge) 
to record quantitatively the complete process of the 
retraction of rubber. However, early work by Exner! (using 
a method by Stefan) is equivalent to a determination of a 
certain average velocity of the tip of a retracting rubber 
sample. Exner measured the velocity of the wave pulse, 
which is simply connected with the velocity of the tip for 
materials with small internal friction like Hevea rubber.’ 
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Fic. 4. Maximum velocity of tip vs. initial elongation for various stocks. 


Simultaneously with our work, but independently, an 
interesting method to measure an average velocity of the 
tip of a sample, over a certain part of the retraction, was 
developed at the Goodyear Research Laboratories.* 
Whereas our method seems to be more suitable for quanti- 
tative work, the Goodyear method seems to be quite 
suitable for factory control purposes of stocks. 

1F, Exner, Ann. d. Physik 153, 62 (1874). 


2 Cf. H. M. James and E. Guth, Phys. Rev. 66, 33 (1944). 
+R. B, Stambaugh, Phys. Rev. 65, 250 (1944). 





Photographic Study of the Retraction of 
Stressed Rubber 


B. ADALBERT Mrowca, C.S.C., S. LEoNARD DART, AND EUGENE GUTH 
Department of Physics, University of Notre Dame, Notre Dame, Indiana 
July 3, 1944 





HE peculiar element by element contraction of rubber, 
discovered by the two-stylus experiments, has been 





Fic. 1. Stroboscopic photographs showing retraction of Hevea gum. 





Fic, 2. Stroboscopic photographs showing retraction of Butyl tread. 
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described.! This behavior was checked by high speed 
photography. Instantaneous photographs of the progress 
of the retraction were taken by a camera in stroboscopic 
illumination. Figure 1 shows a set of such photos for Hevea 
gum. The white marks are equidistant at the initial 
elongation (about 100 percent). The wave pulse and its prog- 
ress are clearly visible. Figure 2 shows a set of photos for 
Butyl tread. Here the pulse is quite diffuse and is reflected 
at the fixed end before the sample has contracted com- 
pletely, leading to an early buckling. Of course, the con- 
traction starts with the tip. Figure 3 shows the differential 
extension, i.e., the extension of each element as a function 
of the distance from the fixed end for two consecutive times 
t and és, respectively. Figure 3 was plotted from the photos 
in Fig. 1 and Fig. 2. For Hevea gum, the progress of the 
well-defined wave pulse is revealed again. For Butyl tread 
the velocity of each point is approximately proportional to 
its distance from the fixed end. 

Instantaneous photos do not give a time scale. Of course, 
the timing may be obtained from our elongation-time 
curves.! Moreover, the timing may be obtained by using 
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Fic. 3. Differential extension as obtained from 
stroboscopic photographs. 


Edgerton’s high speed camera. We have obtained high 
speed movies with timing; the time was recorded on the film 
by a spark at 1/120 second intervals. These movies check 
well the timing obtained by our stylus experiments, namely, 
that the duration of the retraction for about 100 percent 
initial elongation is from 5 to 20 milliseconds according to the 
particular sample used. After the discovery of the element 
by element contraction by our two-stylus experiments, 
instantaneous photos were taken about simultaneously at 
the Goodyear Research Laboratories.” 

A complete theoretical explanation of our findings is 
offered in the following letter by James and Guth.’ Hevea 
gum and Butyl tread represent the two limiting cases of 
small and strong internal friction, respectively. The ex- 
pression for the velocity of the tip: v»=«(E/p)+ agrees® 
rather well with experiment. 


1B. A. Mrowca, S. L. Dart, and E. Guth, Phys. Rev. 66, 30 (1944). 
2R. B. Stambaugh, Phys. Rev. 65, 250 (1944). 
*H. M. James and E. Guth, Phys. Rev. 66, 33 (1944). 
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Theory of the Retraction of Stressed Rubber 


HuBert M. James, Department of Physics, Purdue University, West 
Lafayette, Indiana 


AND 


EuGENeE Gut, Department of Physics, University of Notre Dame, 
Notre Dame, Indiana 


July 3, 1944 


HE striking behavior of stressed rubber in free 
retraction reported by B. Adalbert Mrowca, S. 
Leonard Dart, and Eugene Guth! can be explained by the 
application of the theory of stress-wave propagation ia 
elastic solids. 

We consider a uniform elastic rod of length /, fixed at one 
end, with linear density p in the unstretched state. We 
assume that Hooke’s law holds, with Young’s modulus E, 
and that frictional forces are proportional to the rate of 
deformation. Let x be a measure of distance in the material 
of the unstretched rod, starting at the fixed end, and let u(x) 
be the displacement of the material corresponding to the 
indicated value of x. In the absence of body forces, the 
behavior of the rod in extension and retraction is governed 
by the wave equation 


(Hu /dl) = f(d*u/dx*dt) + E(Pu/dx*). (1) 


If the rod is initially uniformly stretched with fractional 
extension ¢, and is released at time ¢=0, the boundary con- 
ditions are: 
Ou/dx=e, du/dt=0 for all x at t=0, 
u=0 for all ¢ at x=0, (2) 
S(@u/dxdt)+E(du/dx)=0 for x=! when t>0. 


In the limiting case of small friction, f=0, the solution 
becomes, with v= (E/p)3, 
u(x, t) =ex l—x>ut, 
u(x,t)=e(l—vt) l—x<ut<l. (3) 
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According to this solution a pulse of acceleration arises at 
the free end at the moment it is released, and travels, with 
respect to the unstretched material, at a rate v; its rate of 
transmission in space is v, = (e+ 1)v. The velocity of the tip 
of the bar is »; = —du/dt =e. Each portion of the bar is sta- 
tionary until it is reached by this pulse. As the pulse passes 
over it, the stress is completely relaxed, and the material 
reaches its final velocity «. When the acceleration pulse 
reaches the fixed end of the material, after time ¢, =//», the 
bar is unstressed and in uniform motion at the rate ev; the 
original potential energy is completely converted into 
kinetic energy. After this time a solution other than Eq. (3) 
begins to apply. The wave pulse is reflected at the fixed end 
and a pulse of deceleration begins to progress back down 
the bar. In the case of a rod in which transverse motion is 
not prevented, as in the case of a retracting rubber band, 
buckling begins at the moment of reflection of the pulse. 

Deviations from Hooke’s law cause broadening of the 
contraction pulse. Consideration of momentum and energy 
transfer makes it evident at once that an acceleration pulse 
of fixed form and limited extent cannot occur except for a 
Hooke’s law bar. 

In the case of high internal friction the inertia term in 
Eq. (1) may be neglected. The solution then becomes 


u(x, t)Sex exp [—(E/f)t]. (4) 


Equation (4) holds after the wave passed down the sample; 
there is a discontinuity at t=0. 

For treatments of the intermediate case, forced vibra- 
tions, and the effect of a mass loading at the free end, 
reference must be made to a full treatment of the problem 
which is to follow. 


1B. A. Mrowca, S. L. Dart, and E. Guth, Phys. Rev. 66, 30 (1944). 
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Proceedings of the American Physical Society 


MEETING AT ROCHESTER, NEW YORK, JUNE 23-24, 1944 


HE 261st meeting of the American Physical 
Society, being the 1944 Summer Meeting 
in the East, was held in Rochester (New York) 
on Friday and Saturday, June 23 and 24. Our 
meeting was held conjointly with one of the 
American Association of Physics Teachers, and 
comprised the inaugural meeting of our own 
Division of High Polymer Physics. To that As- 
sociation and this Division belongs most of the 
credit for arranging an excellent threefold meet- 
ing, with 238 registrants and 180 in attendance 
at the dinner. The Association had a symposium, 
“Essential physics for a naval career,” and 
another, ‘The training of physicists for in- 
dustry,’”’ in which we were graciously allowed to 
join as co-sponsors; it had also two programmes 
of contributed papers. We had two programmes 
of contributed papers (the abstracts follow) and 
an invited paper, “‘Ultraviolet spectra of biologi- 
cally important molecules,’’ by G. A. Anslow. The 
host institution was the University of Rochester, 
and to Professor F. C. Fairbanks and his depart- 
ment and secretarial staff warm thanks are ex- 
pressed. The weather was merciful for the season. 
The dinner of Friday evening at the Hotel Seneca 
was followed by brief speeches by A. J. Dempster, 
L. McDowell, J. H. Dillon, and L. A. DuBridge, 
and an address, ‘“The new world and the scien- 
tist,’’ by C. E. K. Mees. 

The Division of High Polymer Physics of our 
Society made its debut with an inaugural pro- 
gramme of four scientific sessions comprising 
twenty-one invited papers. The titles of these 
papers are printed below; abstracts and in some 
cases complete texts are to be published elsewhere. 
This programme was arranged by J. H. Dillon, 
M.L. Huggins, W. J. Lyons, and H. F. Mark, who 
with W. F. Busse served as chairmen. A draft of 
by-laws was evolved and adopted by the Division 
in business session and by the Council of the 
Society, and in accordance with its stipulations 
the Council appointed an initial executive com- 
mittee of seven members to serve until the 1944 
Annual Meeting of the Society. These seven are: 
R. B. Barnes, W. F. Busse, P. Debye, J. H. 
Dillon, J. P. Elting, W. J. Lyons, and L. A. 
Wood. The Division of High Polymer Physics is 
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now launched and fully equipped and ready to 
proceed under its own steam. Communications 
should be addressed to its Secretary Dr. W. J. 
Lyons, 2100 Robert E. Lee Boulevard, New 
Orleans 19, Louisiana. 

Other actions by the Council at its meetings of 
Friday and Saturday included the appointment 
of a ‘‘steering committee” of three members (K. 
K. Darrow, A. J. Dempster, R. S. Mulliken) to 
guide the formation of the Division of Spec- 
troscopy authorized at the April meeting, with 
power to add other members to the Committee 
and to arrange symposia; and the appointment 
of a committee of six (S. Dushman, T. A. Read, 
F. Seitz, W. Shockley, S. Siegel, R. Smoluchowski) 
to arrange a symposium and take other steps to 
promote a Division of Physics of the Solid State. 

The Council elected four candidates to Fellow- 
ship (one had previously been a Member) and 
fifty-five candidates to Membership; their names 
are subjoined. 

The Society has lost through death four of its 
Fellows: S. Ballantine (Ballantine Laboratories), 
N. A. Kent (Boston University), Margaret E. 
Maltby (Barnard College), and O. M. Stewart 
(University of Missouri). 

Elected to Fellowship: J. Hillier, H. F. Mark, E. Ott, E. 
M. Purcell. 

Elected to Membership: Bartnoff, Shepard; Barton, David 
Maxwell; Brody, Selma Blazer; Brown, Barremore Beverly; 
Budd, Chester B.; Burns, Preston; Claassen, Howard H.; 
Clark, Orrin H.; Ellis, James; Ferguson, Jasper N.; Gibson, 
Kenneth Alvin; Gilbert, Edgar N.; Graves, J. Dentler; 
Halteman, Eber K.; Harth, Erich M.; Hook, Alonzo Lohr; 
Jastrow, Robert; Kaufmann, Stefan G.; Kessler, Myer; 
Kinnaird, Richard F.; Labin, Edouard; Lanzl, Lawrence 
H.; Levine, Samuel; Lewin, Gerhard; Lopes, Jose Leite; 
Madansky, Leon; Maiden, Robert M.; Mathieson, Ray- 
mond L.; McCarthy, Edward L.; Munier, John Hammond; 
Pardee, Otway O’Meara; Perkin, Richard S.; Pleasonton, 
Frances; Plew, Herman E., Jr.; Prather, Robert W.; 
Proskauer, Eric S.; Rall, Waldo; Rich, Joseph Anthony; 
Schwartz, Don L.; Schwarz, Edward Robinson; Sharer, 
William W.; Shaw, Morton R., Jr.; Sreb, Jules H.; Stephen- 
son, John; Strong, Guy E.; Tallmadge, Francis K.; Thomas, 
Benjamin W.; Vogel, Charles P.; Weltmann, Ruth N.; 
Wetmore, Warren L.; White, Donald W., Jr.; Woods, 
Frank R., Jr.; Wundheiler, Alexander; Urquhart, Robert 
Cameron; Ziomek, Joseph S. 

KarL K. DARROW 
Secretary 
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Symposia of the American Association of Physics Teachers and the American 
Physical Society 
Essential Physics for a Naval Career 

Introduction. COMMANDER F. K. ELper, USN (Ret.), Head of Physics Committee, U.S.1 
Vectors. L1EUTENANT COMMANDER J. A. TIEDEMAN, USNR, Instructor in Physics, U.S.1 
Mechanics. LIEUTENANT COMMANDER L. E. K1InsLer, USNR, Instructor in Physics, U.S.1 
Gyroscopes. LIEUTENANT COMMANDER J. D. Riccin, USNR, Instructor in Physics, U.S.N 
Sound. LIEUTENANT E. R. Pinkston, USNR, Instructor in Physics, U.S.N.A. 

Optics. LIEUTENANT (j.G.) R. A. Goopwin, USNR, Instructor in Physics, U.S.N.A. 


2 zee 
a 


The Training of Physicists for Industry 


Postwar Training of Physicists for Industry. Sau DusHMAN, General Electric Company. 

Training for the Small Industrial Research Laboratory. H. L. Mason, Taylor Instrument 
Companies. 

Training of Physicists for the Optical Industry. Br1an O’Brien, University of Rochester. 

Observations on the Training of Physicists. M. N. States, Central Scientific Company. 

Discussion, opened by J. T. LITTLETON, Corning Glass Company. 


Invited Papers of the Inaugural Programme of the Division of High 
Polymer Physics 


Systems with Superimposed Viscous and Elastic Behavior. A. Topo_sky, Princeton University. 

Transition Phenomena in High Polymers. W. O. BAKER, Bell Telephone Laboratories. 

A Contribution to the Kinetic Theory of Rubber Elasticity. T. ALFrey, Monsanto Chemical 
Company. 

A Recent Development in the Theory of Rubber Elasticity. HusErtT M. JAMEs AND E. Guth, 
University of Notre Dame. 

Electrostatic Properties of Rubber and GR-S. R. S. HAvenniILL, H. C. O'BRIEN, Anp J. J. 
RANKIN, St. Joseph Lead Company. 

Speed of Retraction of Rubber. R. B. SramMBAUGH, MARGARET ROHNER, AND S. D. GEHMAN, 
Goodyear Tire & Rubber Company. 

Rise of Temperature on Fast Stretching of Butyl Rubber. S. L. Dart anp E. Gutn, University 
of Notre Dame. 

Limiting Law of the Reinforcement of Rubber. HuGH M. SMALLWoop, U.S. Rubber Company. 

Relations Between Stress, Strain, and Temperature in a Pure-Gum Vulcanizate of GR-S. 
FRANK L. ROTH AND LAWRENCE A. Woop, National ae of Standards. 

Some Low Temperature Properties of Elastomers. F. S. CONANT AND J. W. Liska, Firestone 
Tire & Rubber Company. 

Physical Properties of Elastomers at Low Temperatures. T. L. LouGHBorovuGn, H. E. Greene, 
AND C. W. Harris, B. F. Goodrich Company. 

Molecular Weights and Weight Distribution of High Polymers Determined from Light 
Scattering. P. M. Dory, B. H. Zimm, anp H. Mark, Polytechnic Institute of Brooklyn. 

Theory of Filler Reinforcement. E. Guth AND HuBert M. James, University of Notre Dame. 

Some Factors Influencing the Brittle Point of High Polymers. R.S. SPENCER AND R. F. Boyer, 
The Dow Chemical Company. 

Kinematic Errors in Bobbin Drives on Roving Frames. S. L. GERHARD, United States Rubber 
Company. 

Control of Elongation in Highly Stretched Tire Cord. H. J. Puitipr anp C. M. Conran, 
Southern Regional Research Laboratory, U. S. Department of Agriculture. 

The Influence of Velocity-Gradient on the Relation Between Viscosity and Concentration 
in Cuprammonium Solutions of Cellulose. W. J. Lyons, Southern Regional Research Labora- 
tory, U. S. Department of Agriculture. 

Orientation of Hydrophile Xerogels in an Electric Field, Part II. S. E. SHepparp anp P. T. 
NEWSOME, Eastman Kodak Company. 

Determination of Molecular Weights of High Polymers by Solubility Methods. D. R. Morry 
AND J. W. TAMBLYN, Eastman Kodak Company. 

Comparison of the Structures of Stretched Linear Polymers. M. L. HuGGins, Eastman Kodak 
Company. 

Fractionation, Viscosity, and Osmotic Pressure Studies on Cellulose Esters. J. W. TAMBLYN, 
R. L. TicHener, D. R. Morey, AND R. H. WAGNER, Eastman Kodak Company. 











36 AMERICAN PHYSICAL SOCIETY 


Abstracts of Contributed Papers 


Al. X-Ray Coincidences From Te'*! (125 Days).* Rosa- 
LYN S. YALOW AND M. GOLDHABER, University of Illinois.— 
Te! (125d) which decays by K-electron capture into Sb! 
can be produced by deuteron bombardment of antimony- 
The radiations emitted by Te are very complex, consist. 
ing of x-rays, internal conversion electrons, and y-rays. By 
critical absorption and fluorescence studies, O’Neal and 
Scharff-Goldhaber found not only the Ka lines of Sb but 
also x-rays of somewhat shorter wave-length (‘‘tellurium 
like’). They interpreted these x-rays as due to a transition 
in the Sb atom after both K electrons had been removed, 
one by K-electron capture and the other by a succeeding 
internal conversion of a y-ray. We have obtained further 
evidence strengthening this interpretation by observing 
Geiger counter coincidences between the x-rays emitted 
during the decay of Te™. The X-X coincidence rate 
observed was about 10-* per x-ray counted in a single 
counter. When the x-ray rate in a single counter is reduced 
to about half by interposing different absorbers of suitable 
thickness between the Te sample and the counters, it is 
found that Cd foils, which absorb preferentially the shorter 
wave-length member of the x-ray pair, reduce the coinci- 
dence rate more strongly than either Ag or Al foils. From 
the relative intensity of the “double hole’’ x-ray transition 
to that of the “single hole’”’ x-ray transition, one may 
obtain an estimate of the absolute mean lifetime of inter- 
nally converted nuclear excited states by using for com- 
parison the mean lifetime for which a hole in the K-electron 
shell exists. This time can be calculated from the known 
widths of x-ray lines. One may conclude that the excited 
state of the Sb™ nucleus which gives rise to strong internal 
conversion lines has a mean lifetime +2 X10~"’ sec. 


* Read by title. 


A2. Slowing Down of R Neutrons into C Neutrons in 
Lead. C. O. MUEHLHAUSE AND M. GOLDHABER, Uni- 
versity of Illinots—When fast neutrons are slowed down 
inside paraffin they pass through an energy region above 
the Cd absorption limit where they may be referred to as 
R neutrons. After a few further collisions with protons they 
become C neutrons. It is of interest to study with what 
probability R neutrons can be slowed down into C neutrons 
by interaction with heavy elements. A previous attempt to 
detect this transformation was unsuccessful. With a more 
sensitive arrangement we have now been able to detect 
such an effect in lead and other elements. Its magnitude 
depends to some extent on geometrical and other factors. 
With one arrangement the slowing down effect from R into 
C neutrons per lead atom was found to be 5 percent of that 
observed per hydrogen atom. It is unlikely that more than 
a small part of the effect is due to hydrogen contamination 
in the lead. The effect appears to be too large, however, to 
be accounted for by excitation of lattice vibrations alone. 
No effect was observed with our arrangement with either 
Ra-a-Be or Ra-y-Be sources when no paraffin was present 
to slow down the primary neutrons into the R region. 


A3. Nuclear Excitation of Cadmium. Marce..us L, 
WIEDENBECK, University of Notre Dame.—By x-ray irradia- 
tion of self-quenching argon-ether counters having cad- 
mium cathodes, the excitation curve for this element was 
obtained up to 2.9 mev. The half-life+? was found to be 
48.7+0.3 minutes. The threshold for the excitation is 1.25 
+0.03 Mev and activation levels were found at 1.68 Mev, 
2.08 Mev, and 2.56 Mev. The electron excitation curve 
shows the same threshold and activation levels which were 
found by x-ray excitation. Absorption of the conversion 
electrons by aluminum foils indicates that the energy of the 
metastable state is 195 kv. 


1 Dode and Pontecorvo, Comptes rendus 207, 287 (1938). 
2 J. R. Feldmeier and G. B. Collins, Phys. Rev. 59, 937 (1941). 


A4. The Electron Affinity of Bromine and a Study of 
Its Decomposition on Hot Tungsten. P. M. Doty anp 
J. E. MAvER, Columbia University.—The electron affinity 
of bromine has been measured by determining the ratio of 
ions to electrons leaving a hot tungsten surface in contact 
with bromine gas at low pressure. The value obtained is 
80.5+0.4 kcal./mole. The measurements, which extend 
over a range of 700° in temperature and 80-fold in pressure, 
show no trend in the electron affinity and demonstrate 
that the decomposition of bromine under these conditions 
is a first-order reaction. The assumption of unit accommo- 
dation coefficient is shown to be correct. The agreement 
between the value here obtained and the value calculated 
from energy cycles substantiates the improved method of 
calculating lattice energies. 


AS. Analysis of the Raman Spectrum of a New C;H, 
Compound. Forrest F. CLEVELAND AND M. J. Murray, 
Illinois Institute of Technology—Murray and Stevenson! 
recently reported Raman frequencies and relative intensi- 
ties for a new CsHs compound and presented five possible 
structures for the molecule. The method of preparation as 
well as the values of the C—H stretching frequencies indi- 
cated that the compound was spiropentane. In the present 
investigation, the depolarization factors of the Raman lines 
have been obtained and the expected number of polarized 
and depolarized frequencies, both for the entire molecule 
and for the carbon skeleton (Av less than 1050 cm), have 
been worked out for each of the five structures by use of the 
group theory formulas. If it is assumed that all lines in- 
volved in accidental degeneracies have been resolved, the 
structure giving the best correlation with the experimental 
results is one having the symmetry Dy. However, if un- 
resolved accidental degeneracies do exist, as seems quite 
probable for the spiropentane structure since there are a 
number of similar groups in the molecule, the spiropentane 
structure of symmetry C2, would give the best correlation. 
A calculation of the frequencies would aid in deciding be- 
tween the two structures. The large number of approximate 
coincidences between the Raman and infra-red frequencies 
is additional evidence in favor of the spiropentane structure. 


! Murray and Stevenson, J. Am. Chem. Soc. 66, 314 (1944). 
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A6. The Formation of a Zone Plate by Newton’s Rings. 
RoBERT AyRES Woopson, Eastman Kodak Company.— 
The approximate equation for the radius y’ of the Nth ring 
of a zone plate is: y’*= Nfx= Ny,"2, where f =focal length, 
\= wave-length, y:’=radius of first ring. The approximate 
equation for the radius y”’ of the Nth Newton ring is: 
y’""*= Nd/n(Ci— C2) = Ny", where } = wave-length, m =in- 
dex of refraction of thin film, C;=curvature of one surface 
of thin film, C2=curvature of other surface of thin film, 
y:” =radius of first ring. Since the radius of the Nth ring is 
proportional to ¥ N for both the zone plate and Newton's 
rings, it occurred to the author to use Newton's rings for 
the formation of a zone plate. On combining the above 
equations, we obtain for the power P of such a zone plate, 
P=n(C,—C:), or for an air film between the lens surfaces, 
P=(C,—C;. For the combination of a plano-parallel plate 
and the convex surface of a zero-power meniscus lens, the 
focal length is equal to the radius of the convex surface. 


A7. Some Basic Thermodynamic Considerations. J. L. 
Finck, College of the City of New York.—In previous papers 
the writer has endeavored to show that from a phenomeno- 
logical point of view there is justification for including 
metastable states within classic thermodynamics, particu- 
larly if catalytic phenomena are considered. In doing so, a 
much broader view may be taken of the thermodynamic 
behavior of systems. However, criteria based exclusively on 
the first two laws of thermodynamics define equilibrium in 
such manner that metastable states are excluded. In the 
present paper a critical examination is made of develop- 
ment of the second law, and a fallacy is found wherein the 
metastable states are excluded from the very start—in the 
treatment of Carnot’s cycle. Logical approach requires 
that we establish the independent variables in advance, and 
it is found that in general an equation of state should 
include more than two independent variables. Such an 
equation will include metastable, as well as ordinary equi- 
librium states, and will apply to systems consisting of one 
or several phases. In enlarged systems it is shown that 
entropy is no longer an integrable function, and therefore 
cannot play the role it does in classic thermodynamics. It is 
suggested how an enlarged equation of state may explain 
thermodynamically the catalytic behavior of systems where 
intermediate products are formed. 


Bl. Knudsen Type Vacuum Gauge in Technical High 
Vacuum Work. A. F. TuRNER AND O. A. ULLRICH, JR., 
Bausch & Lomb Optical Company.—The Knudsen vacuum 
gauge has not found widespread application. Apparently 
the principal objection hitherto has been the need for a 
vibration-free support. This restriction can be materially 
relaxed in the design of DuMond and Pickels' by the addi- 
tion of a small auxiliary internal damping vane moving in 
a horizontal plane between permanent magnets. External 
permanent magnets damp the impact vane. Gauges so 
constructed and mounted rigidly on the base plates of large 
bell jar systems have given satisfactory performance in 
evaporation work under adverse floor vibration conditions. 
The gauge is practically linear through the 10-* mm Hg 
range using a sensitivity of 1X 10-* mm Hg per cm deflec- 
tion at 23-cm scale distance. Operation in-higher and lower 
ranges is obtained by changing the heater temperature. 


The gauge is not harmed if accidentally turned on at 
atmospheric pressure—a definite advantage in a technical 
instrument. The performance of the gauge during the eva- 
poration of some metals and salts will be described. 


assy W. M. DuMond and W. M. Pickels, Jr., Rev. Sci. Inst. 6, 362 


B2. Solution of the Diffusion Equation Applicable to the 
Edgewise Growth of Pearlite. W. H. Branpt, Westing- 
house Electric & Manufacturing Company.—The diffusion 
equation is transformed to a set of coordinates moving with 
the pearlite interface and a solution applicable to the 
problem obtained in the form of an infinite series of terms. 
Using the first three terms, the edgewise velocity of pearlite 
growth is calculated for a plain carbon eutectoid steel using 
data most of which are obtained by extrapolation. The 
values obtained show reasoénable agreement with values 
for the rate of pearlite nodule growth determined by Hull, 
Colton, and Mehl.* The velocity increases with decreasing 
temperature, as expected, and it is shown that this is due 
to the change in the solubility of ferrite and cementite in 
austenite with temperature. The theory predicts curved 
ferrite and cementite interfaces and the carbon concentra- 
tion in austenite is shown to vary across each of these 
interfaces. : 


* F.C. Hull, R. A. Colton, and R. F. Mehl, Trans. A.I.M.E. 150, 185 
(1942). 


B3. The Mechanics of the Metal Cutting Process. II. 


“Oblique Straight-Edged Tool and Continuous Chip. M. 


EUGENE MERCHANT, Cincinnati Milling Machine Company. 
—The principles employed in the derivation of the me- 
chanics of the simple case previously discussed! have been 
extended to an analysis of the geometry of chip formation 
and of the force system occurring in the more general case 
of a straight-edged tool with cutting edge oblique to the 
direction of relative motion of tool and workpiece. It is 
found that the direction of chip flow and the direction of 
shear are no longer perpendicular to the cutting edge as 
they are in the simple case. As a result the force system is 
three-dimensional. However, it is again possible to obtain 
simple relationships between the readily measurable force 
components and geometrical variables, and such basic 
quantities as strain in the chip, friction between chip and 
tool, shear strength of the metal, etc., thus making possible 
the analysis of actual machining operations in terms of 
these basic quantities. 


1M. E. Merchant, Bull. Am. Phys. Soc. 19, No. 2, 7 (1944). 


B4. Experimental Investigations of the Mechanism of 
Friction. B. W. SAKMANN, Massachusetts Institute of Tech- 
nology.—In all friction experiments material is exchanged 
between the sliding surfaces even under small loads of a 
few grams. This fact seems to indicate that such an ex- 
change of matter might be of fundamental importance in 
all frictional phenomena. The development of a radioactive 
tracer method made it possible to investigate this transfer 
of matter. Spherical riders were slid over an activated base 
surface and were tested for the presence of radioactive 
material. The amount of material transferred depends upon 
various parameters. The amount of base metal adhering to 
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the rider was found to be proportional to the normal load 
and the distance traveled. The transferred material in- 
creased with the surface roughness if the rider was harder 
than the base. If, on the other hand, the rider was softer 
than the base, surface finish was of minor importance. For 
steel samples the deposited material was inversely propor- 
tional to the Brinell hardness. For riders of different metals 
of the same hardness and surface finish the amount of 
transferred material increased with the solid solubility of 
the base metal in that of the rider. The exchange of metal 
was reduced by lubrication. The reduction depended on 
the normal load and, for the same load, on the material of 
the rider. 


BS. Tensile Strength of Liquids. ALLEN KiNG, Dart- 
mouth College—When a tension is applied slowly to a liquid 
column, rupture occurs with the formation of a cavity. 
The critical tension for rupture probably fluctuates about 
a mean. The smallest cavity in which vapor can exist in 
equilibrium within a liquid at absolute temperature T may 
be assumed to have a diameter equal to the average dis- 
tance between molecules of vapor at the vapor pressure for 
this temperature. This cavity is in unstable equilibrium. If 
q represents the critical tension (tensile strength) in at- 
mospheres, then g = 14.67(~/T)!, where y is the free surface 
energy in ergs per cm? and is the vapor pressure in cm of 
Hg. At 20°C g has a value of 194 atmos. for water, 80 atmos. 
for ethyl alcohol, and 132 atmos. for diethyl ether. Accord- 
ing to this relation gq increases with temperature to a 
maximum value. 


B6. Ponderomotive Forces on Matter in Electric and 
Magnetic Fields. Fer1x EnxreENunAFT, New York City.—In 
a recent Letter to the Editor of The Physical Review it was 
stated that beside the force of gravity and the magnetic 
action of electric currents (Oersted, Biot-Savart, H. A. 
Lorentz), f.=D+1/c[v., H], there is a third general force, 
the electric action of magnetic currents, f, = B+1/c[vm, E). 
However, if one adds to the electrodynamic equations the 
term of the magnetic current, following the conception of 
Stefan and Heaviside, it seems that the forces f, and fn 
cannot explain all the observed movements. If there is no 
angle between the direction of the electric current and the 
magnetic field, and no angle between the direction of the 
magnetic current and the electric field, the second terms of 
the above equations become nil. But it has been observed 
that, when the homogeneous electric and magnetic fields 
are parallel and in action simultaneously, there can be a 
circular movement of matter. One suspects that some of the 
phenomena reported by S. P. Thompson and others are of 
the same type, although the homogeneity of the fields they 
employed was not considered by them. Experiments in 
gases as well as in liquids will be described where these 
movements are to be seen in homogeneous fields. Con- 
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clusions can only be drawn under the preliminary assump- 
tion that there is no interdependency between H and E 
other than that expressed in the equations mentioned 
above. 


Read by Title 


Tl. The Electrical Resistivity of Alpha-Phase Copper- 
Zinc Alloys at Low Temperatures. H. A. FarrBANK, Yale 
University.—A systematic investigation of the resistivity 
(e) of the copper-zinc alpha-phase solid solution system 
has been made in the temperature range 14.3°K to room 
temperature. Measurements were taken on two sets of one 
millimeter diameter wires. The first was annealed at 700°C 
and quenched; the second, cold-drawn through a 75 percent 
reduction and left in the hard-drawn state. For the an- 
nealed set, the resistivity varied linearly with the tempera- 
ture above 70°K and with a higher power of T at low 
temperatures, approaching a constant value below 20°K. 
Above 70°K, dp/dT was found to increase linearly with 
atomic percent solute in variance with Matthiessen’s rule. 
The residual resistivity increased with percent added zinc, 
but not in agreement with Nordheim’s theory,' unless it 
be assumed that the effective number of free electrons per 
atom increases with added solute atoms or that the speci- 
mens were not true random solid solutions. The hard-drawn 
specimens showed properties similar to the annealed group 
except for an increase in the resistivity due to residual 
strains. 


1L. Nordheim, Ann. d. Physik 9, 664 (1931). 


T2. The Electrical Resistivity of Alpha-Phase Copper- 
Tin Alloys at Low Temperatures. H. A. FAIRBANK, Yale 
University.—An investigation similar to that described in 
the previous abstract for the Cu-Zn alpha-phase alloys 
has been carried out for the Cu-Sn alpha-phase alloy system 
in the temperature range 14.3°K to room temperature. The 
dependence of the electrical resistivity on temperature and 
on percent added tin was similar to that found for the 
Cu-Zn system. However, the cold work produced some 
unexpected effects. Although the expected increase of re- 
sistivity with residual strain occurred in the higher per- 
centage tin specimens, the wires of zero to three atomic 
percent tin, which were hard-drawn through a 75 percent 
reduction and left in the strained state, showed a lower 
resistivity at low temperatures than did the corresponding 
annealed, quenched specimens. This anomaly increased as 
the temperature was lowered. This unusual behavior may 
be due to order in the Cu-Sn specimens of about two atomic 
percent tin, induced either during the cold working process 
or during the air cooling which followed the anneal at 
750°C prior to the 75 percent reduction. It is also possible 
that a phase separation may have occurred as a result of 
the extensive cold work. Further investigation appears 
highly desirable. 
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